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1.  GERESS  STATUS  REPORT 
December  1993  -  May  1994 


Michael  L.  Jost 


1.1  General 

The  German  Experimental  Seismic  System  (GERESS)  has  been  operated 
by  Ruhr-University  Bochum.  The  GERESS  array  is  located  in  southeastern 
Germany  (Bavarian  Forest)  near  the  brader  to  Austria  and  the  Czech  Republic. 
The  array  consists  of  25  stations  with  vertical-component  short-period  Teledyne 
Geotech  GS-13  instruments  sampled  at  40  Hz.  Four  of  the  sites  include  hor¬ 
izontal  component  instruments.  At  the  key  station  of  the  array,  GEC2 
(48.84511  N,  13.70156  E,  1132  m),  a  STS-2  is  sampled  at  80  Hz  (broad-band 
element).  The  geometry  of  the  array  is  based  on  concentric  rings  providing  an 
overall  aperture  of  about  4  km  (Harjes,  1990).  GERESS  was  installed  in  1990 
and  the  number  of  recording  elements  increased  gradually  over  that  year.  The 
array  became  fully  operational  in  January  1991.  Data  from  the  array  have  been 
continuously  transmitted  to  the  NORSAR  data  center  in  Kjeller  and  to  the  Insti¬ 
tute  of  Geophysics,  Ruhr-University  Bochum  via  64-kbit  lines. 

This  repon  summarizes  the  technical  performance  of  GERESS  and  the 
Bochum  data  center,  and  continues  similar  descriptions  given  previously  (lost, 
1991a;  1992a;  1993a,b).  This  summary  is  based  on  monthly  status  reports  that 
are  available  upon  request. 


1.2  GERESS  Data  Center  at  Ruhr-University  Bochum 

The  Institute  of  Geophysics  of  Ruhr-University  Bochum  has  operated  an 
experimental  on-line  processing  system  for  GERESS  data  (Jost  1992b;  Harjes  et 
al.,  1993).  This  system  uses  software  developed  at  NORSAR  data  center  (Myk- 
keltveit  and  Bungum,  1984;  Fyen,  1987).  The  automatic  analysis  has  worked 
very  stable  and  results  are  available  for  all  GERESS  data  in  this  reporting 
period. 

The  on-line  processing  consists  of  3  steps:  detection,  fk-analysis,  and  loca¬ 
tion.  The  first  stage  of  the  on-line  processing  accesses  data  in  30  second  seg¬ 
ments  and  runs  a  STA/LTA  detector.  The  detector  presently  recognizes  an 
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onset  if  the  STA/LTA  ratio  for  a  filtered  trial-beam  exceeds  a  given  threshold 
ranging  between  3.5-4.0,  Since  September  1993,  horizontal  beams  have  been 
included  to  improve  on  S-type  detections.  The  next  step  of  the  on-line  process¬ 
ing  is  the  tran^oimation  of  a  filtered  data  segment  at  each  onset  time  (derived 
from  the  detection  time)  into  the  frequency-wavenumber  domain.  As  a  result, 
the  slowness  and  back-azimuth  of  the  phase  is  determined.  From  the  slowness 
infftmnaiion,  seismic  phases  are  identified.  The  final  step  of  the  on-line  data 
processing  is  the  location  of  events.  The  seismic  phases  as  identified  in  the  fk- 
analysis  are  associated  to  events  in  this  step.  From  the  arrival  time  difference 
of  regional  phases,  the  distance  to  the  epicenter  is  determined  from  the  Jeffreys 
-  BuUen  travel  time  tables  for  regional  seismic  phases.  Together  with  a  mean 
back-azimuth,  the  epicenter  locations  of  local  and  regional  events  are  deter¬ 
mined.  Correcting  for  outages,  20  events  have  been  automatically  located  each 
day  on  average  (Table  1-1).  Figure  1-1  shows  the  automatic  event  locations  for 
Dec.  1993  -  May  1994.  Local  magnitudes  aut(»natically  determined  at  Bochum 
(e.  g.  lost,  1993b)  i^ipeared  low  by  at  least  0.5  magnitude  units  compared  to 
values  published  by  other  European  observatories.  A  new  amplitude-distance 
curve  has  been  implemented  on  March  1,  1994.  A  detailed  description  will  be 
given  elsewhere. 

Alert  messages  for  strong  teleseismic  and  regional  events  have  been 
automatically  sent  to  NEIC  (Golden,  CO)  and  other  interested  institutions  in 
near  real-time.  44  tunes,  a  GERESS  alert  message  for  a  teleseismic  onset  has 
been  used  in  the  NEIC  alerts  in  the  lepmting  period;  and  19  times  by  the 
EMSC.  For  regional  events,  a  fast  earthquake  information  system  (FEIS)  has 
been  developed  (Schweitzer,  Schulte-Theis,  Jost)  and  operated  at  Bochunt 
Starting  from  a  GERESS  automatic  detection  above  a  certain  threshold,  the  rou¬ 
tine  automatically  accesses  the  detection  logs  of  the  German  Regional  Seismic 
Network  (GRSN).  Using  all  available  detections  from  this  network,  the  routine 
proceeds  to  calculate  a  new  location  which  is  subsequently  e-mailed  to 
interested  institutions.  The  time  difference  between  origin  time  and  transmis¬ 
sion  of  an  alert  message  was  42  minutes  for  more  than  84  %  of  all  FEIS  alerts 
(Jul  93  -  Oct  93).  The  shortest  time  was  24  minutes,  the  longest  1  hour.  This 
rime  delay  was  mainly  detetinined  by  the  availability  of  communication  lines. 

The  specific  project  of  re-locating  regional  events  by  a  human  analyst  was 
limited  to  1993.  From  January  1994,  analysts  have  routinely  reviewed  regional 
first  arrivals  together  with  teleseismic  events.  For  this  reporting  period,  parame¬ 
ter  data  (onset  time,  amplitude,  period,  azimuth,  and  slowness)  for  9783  phases 
(i.e.,  54  phases  per  day)  have  been  sent  to  NEIC  (Golden,  C!0),  ISC,  EMSC, 
and  other  institutions. 

GERESS  is  an  open  station  and  results  of  the  on-line  processing  have 
been  sent  to  interested  institutions  via  e-mail  (e.g.  Bundesanstalt  fur 
Geowissenschafren  und  Rohstoffe  (BGR),  EMSC,  Frankfurt  University, 
Geoforschungszentrum  Potsdam  (GFZ),  Grafenberg  Array,  ING,  KTB  array, 
LDG,  Leipzig  University,  Stuttgart  Uruversity,  NORSAR,  Oklahoma  Gcol.  Sur¬ 
vey,  ORFEUS,  and  seismological  institutes  in  Austria,  CZ,  Hungary,  Slovenija, 
Switzerland;  YKA).  In  addition,  the  On-line  Processing  Display  Manager  (Jost, 
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1991b)  has  been  upgraded  several  times  to  make  the  diverse  local  data-bases 
available  for  near  r^  time  information  on  parameter  data  (Fig.  1-2).  The  Data 
Request  Manager  (DRM)  by  Krake  Inc.  (under  contract  from  BGR,  Hannover) 
has  been  operated  in  trial  mode.  Its  performance  has  remained  insufficient  and 
a  frequently  announced  upgrade  has  been  waited  for.  Eventually,  the  DRM 
should  enable  the  transmission  of  data  (GSE,  SEED)  and  interface  to  the  Ger¬ 
man  Regional  Network  of  Broad-Band  Stations. 


Table  1-1 

Statistics  on  the  GERESS  Bochum  On-line  Processing  for  Dec  1993-May  1994 

number  number/day 


detections 

63608 

349 

f-k  analyses 

50531 

278 

locations 

3690 

20 

1.3  GERESS  Array  Hub 

The  GERESS  array  has  shown  stable  operation  from  Dec.  1993  -  May 
1994.  Table  1-2  shows  the  uptime  of  the  system  including  the  Bochum  data 
center  (uptimes  directly  correspond  to  the  data  archived  in  Bochum).  The  aver¬ 
age  uptime  was  98.9  %,  i.e.,  the  total  outage  amounts  to  1  day  23  hours.  The 
highest  uptime  (99.99  %)  was  observed  in  March  1994;  the  lowest  uptime  (95.3 
%)  occurred  in  Dec.  1993  due  to  the  sluggish  response  of  the  German  Telekom 
in  replacing  a  data  noodern  damaged  by  thunderstorm.  Table  1-3  further 
specifies  the  causes  of  the  observed  outages.  The  technical  status  of  GERESS 
is  summarized  in  chronological  order  in  Appendix  1-1.  Results  of  a  recent 
instrument  calibration  are  given  in  Appendix  1-2. 

Table  1-2 

GERESS  uptime  Dec.  1993  -  May  1994  (incl.  Bochum  data  center) 


Month  uptime 
% 


Dec 

95.30 

Jan 

99.95 

Feb 

99.96 

Mar 

99.99 

Apr 

99.96 

May 

98.42 
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Table  1-3 

Causes  of  GERESS  downtimes  Dec.  1993  -  May  1994  (incl.  Bochum  data  center) 
(The  total  downtime  was  171922  sec  =  Id  23h  4Sm  22  sec) 


cause 

downtime 

downtime 

[sec] 

[%] 

Bochum  data  acq.  software 

17075 

9.9 

Bochum  data  acq.  software  (maintenance) 

0 

0 

Bochum  data  acq.  workstation 

5 

0 

Bochum  CIM 

0 

0 

Bochum  power  outage 

41726 

24.3 

total  outage  Bochum 

58806 

34.2 

Communication  line 

75685 

HUB  lAC  crashes  (1) 

98 

0.1 

HUB  thunderstorm  damage 

19774 

11.5 

HUB  maintenance 

3659 

2.1 

HUB  power  outage 

13769 

8.0 

total  outage  HUB 

37300 

21.7 

other  causes 

131 

0.1 

GERESS  has  used  data  acquisition  equipment  from  Teledyne-Geotech  (i.e. 
RDAS-200,  lAC)  which  is  no  longer  in  production.  GERESS  is  a  prototype 
installation  and  inherent  problems  have  remained  (e.g..  Golden  et  al,  1991). 
Due  to  the  severity  of  these  problems,  we  will  continue  to  list  them  to  remind 
the  community  of  possible  problems  inherent  in  GERESS  data. 

-  The  RDAS-2(X)  occasionally  sends  electronic  noise  instead  of  seismic  data 

after  power  up  or  self-reset  This  state  can  be  cleared  only  by  manual 
intervention  (reset). 

-  If  the  Ipps  timing-signal  is  bad  (optical  modems,  fiber  optic  link),  or  the 

clock-board  in  the  RDAS-200  deteriorates  before  completely  failing,  or 
after  any  self-reset,  or  after  any  power  cycle,  data  can  be  desynchronized 
by  one  second  or  by  less  than  a  second  (mostly  0.2  seconds).  Only  an 
array  installation  of  RDAS-200’s  has  the  chance  to  detect  this.  Each  day,  a 
careful  inspection  of  GERESS  records  has  been  performed  to  insure  proper 
synchronization.  After  a  desynchronization  has  been  detected,  one  or 
several  manual  power  cycles  usually  synchronizes  the  station. 

-  The  RDAS-200  is  very  sensitive  to  spikes  in  input  voltage  (thunderstorms) 

causing  outages  or  self-resets. 
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-  The  RDAS-200  occasionally  sends  data  with  crarupted  time  stamps  in  the 

future. 

-  Meteorological  data  are  not  available  for  GERESS. 

-  GERESS  runs  cm  DCF-77. 

At  BcKhum  data  center,  we  have  continued  to  keep  GERESS  operational  to 
provide  the  community  with  useful  data.  We  hope  to  gather  sufficient  spare 
parts  to  keep  the  current  hardware  running  as  long  as  practical.  For  general 
scientific  purposes,  this  situation  appears  acceptable.  But  for  the  special 
requirements  in  a  treaty  situation,  we  recommend  an  upgrade  of  the  present 
hardware  to  a  more  robust  system. 

1.4  Conclusion 

GERESS  is  the  most  sensitive  seisnx>logical  station  in  Central  Europe  for 
regional  events.  The  teleseismic  detecticHi  threshold  of  GERESS  ranges  glo¬ 
bally  between  m),  4.0  -  4.5  (this  repmt). 

At  the  GERESS  data  center  B<x:hum,  the  observatory  routine  has  included 
data  archival  and  automatic  on-line  analysis.  Ftxxn  Dec.  1993  to  May  1994,  20 
regional  events  have  been  automatically  located  each  day  on  average.  The 
GERESS  data  center  Bochum  is  the  only  German  station  that  provides  NEIC 
and  interested  institutions  with  alert  messages  (teleseismic  and  regional). 

In  Bochum,  extensive  parameter  data  for  54  regional  and  teleseismic 
phases  have  been  interactively  reviewed  each  day  on  average.  Results  of  this 
re-analysis  of  GERESS  data  have  been  used  by  NEIC,  ISC,  EMSC,  and  other 
interested  institutions.  The  Bochum  data  center  is  an  open  station  similar  to  the 
German  Regional  Network  of  Broad-Band  Stations  (GRSN).  In  addition, 
Bochum  serves  as  data  center  where  various  bulletins  (e.g.,  alerts  from  NEIC, 
EMSC,  NORSAR,  SED,  and  ING;  onset  lists  from  GRSN,  Austrian  stations, 
Czech  stations,  and  YKA;  locations  from  NEIC  (QED,  FDE  (weekly),  PDE 
(monthly)),  IDC  (GSETT-3),  IGN,  IMS,  LDG,  and  SED;  fast  CMT’s  from  Har¬ 
vard,  USGS,  and  Tokyo  University)  are  available  for  the  interested  user. 

The  GERESS  data  center  Bochum  has  aided  the  CSS  to  test  the  experi¬ 
mental  GSETT-3  International  Data  Center  (IDC  VO,  since  May  1993).  Within 
this  concept,  GERESS  served  as  a-station  and  acted  as  a  node  to  access  data 
from  the  GRSN,  serving  as  ^stations. 

The  GERESS  data  center  Bochum  has  been  responsible  for  GERESS 
maintenance.  The  average  GERESS  uptime  which  includes  the  Bochum  data 
center  has  increased  to  99.65  %  during  the  first  5  months  of  1994.  For  com¬ 
parison,  GERESS  uptimes  were  97.2  %  in  1993,  96.9  %  in  1992,  and  91.5  % 
in  1991.  This  again  demonstrates  the  high  level  of  reliability  of  the  installation, 
its  successful  maintenance,  and  the  high  motivation  of  key  personnel. 
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Appendix  1-1:  Chronological  Status  List  GERESS 


Dec.  1993: 

-  outage  HUB  (Telekom  data  modem  damaged  by  thunderstorm) 

-  power  outage  HUB  (power  line  damaged  during  storm) 

-  power  outage  on  UPS  university  (fuse  blown  by  short  elsewhere) 

-  2  min  gap  possibly  due  to  power  outage  and  defective  USV  (HUB) 

-  Dec.  13  -  14,  Dm.  27  -  Jan  1,  maintenance  visits  at  HUB. 

-  data  link  GERESS  -  Kjellcr  and  GERESS  -  Bochum  switched  from  land-lines 
(German  Telekom)  to  satellite  (Spaceline  Inc.). 

Jan.  1994: 

-  2  min  dafa  gap  possibly  due  to  power  outage  and  defective  USV  (HUB) 

-  2  thunderstorms  at  HUB  (3  times  2  min  outages  array,  station  outages) 

Feb.  1994: 

-  temporary  malfunction  of  STS-2  broad-band  instrument 
Mar.  1994: 

-  2  min  data  gaps  due  to  power  outage  and  defective  USV  (HUB) 

-  implementation  of  new  amplitude-distance  curve  for  magnitude  calculation 

-  modification  of  system  kernel  m  data  acquisition  workstation 

Apr  1994: 

-  2  min  data  gaps  due  to  power  outage  and  defective  USV  (HUB) 

-  defective  preamplifier  C2  (broad-band  element) 

-  5  thunderstorms  at  HUB  (station  outages) 

-  Apr.  10-14  maintenance  visit  at  HUB. 

-  implementation  of  German  3.  level  seismotectonic  names  of  lASPEI  regionaliza¬ 
tion  (Leydecker  1994) 

May.  1994: 

-  3  thunderstorms  at  HUB  (array  and  station  outages) 

-  power  outage  HUB  (maintenance  by  electricity  company) 

-  May  2-6  maintenance  visit  at  HUB  (preamplifier  C2  exchanged,  calibration 
problems  fixed  on  B3,  C7,  D6,  D7,  D8,  D9,  replacement  of  defective  USV) 

-  spontaneous  desynchronization  on  D7 

-  temporary  malfiuiction  of  STS-2  broad-band  instrument 

-  new  USV  defective  (manual  intervention  required  after  surge  or  extended  power 
outage) 
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Appendix  1-2:  GERESS  Instrument  Calibration 


Each  day,  the  GS-13’s  of  GERESS  are  calibrated  using  a  1  Hz  sinusoidal  input  voltage 
(GS-13  4500  mV)  fed  to  the  calibradcm  coil. 


cal.  tiine 

channel 

K<m 

3:00:50 

AO.  A-ring  (verticals).  C2-sz 

300- 

3:01:00 

D6 

3:01  - 

3:01:50 

A2-sn 

3:02- 

3:02:50 

A2-ae.  B  -  ring 

3K)4- 

3K)4:50 

C1.C3 

3K)6- 

3:06:50 

Dl-sz.  D2.  D3 

3K)7- 

3:07:50 

Dl-sn.  D4-sn.  D7-sa 

3.-08  - 

3:08:50 

Dl-se.  D4-SC.  D7-sc.  C4.  C5.  C6.  C7 

3:10- 

3:10:50 

D4-SZ.DS 

3:12- 

3:12:50 

D7-8Z.  D8.  D9 

Calibration  on  May  26,  1994  (for  D6  from  May  8;  318273  ents  are  expected) 


channel 

ent 

deviation  in  % 

GEAO  sz 

317123 

-0.361 

GEArsz 

319254 

0.308 

GEA2  sz 

311215 

-2.22 

GEA3  sz 

316905 

-0.430 

GEBl  sz 

319442 

0.367 

GEB2  sz 

316260 

-0.632 

GEB3  sz 

380812 

19.65 

GEB4  sz 

320425 

0.676 

GEB5  sz 

308851 

-2.96 

GECrsz 

314933 

-1.05 

GEC2  sz 

312405 

-1.84 

GEC3  sz 

315325 

-0.926 

GEC4“sz 

311943 

-1.99 

GEC5  sz 

312094 

-1.94 

GEC6  sz 

315217 

-0.960 

GEC7  sz 

309905 

-2.63 

GEDl  sz 

297258 

-6.60 

GED2  sz 

324986 

2.11 

GED3“sz 

318497 

0.070 

GED4  sz 

300922 

-5.45 

GED6  sz 

317980 

-0.09 

GED5'‘sz 

312444 

-1.83 

GED7"sz 

289007 

-9.20 

GED8  sz 

316283 

■0.625 

GED9“sz 

320454 

0.685 

fa 

318844 

0.179 

GEA2  se 

315890 

-0.749 

GEDl  sn 

324707 

2.02 

GEOrse 

314142 

-1.30 

GED4“sn 

312861 

-1.70 

GED4“se 

313429 

-1.52 

GED7  sn 

309456 

-2.77 

GED7“se 

308842 

-2.96 

8 


GERESS  BOCHUM:  AUTOMATIC  LOCATIONS  DEC.  1  993  -  MAY  1994 


Figure  1-1: 
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Figure  1-2:  The  basic  menu  of  the  On-line  Processing  Display  Manager  at  the 
GERESS  data  center  Bochum 


******************** **************************************************** 


★ 

★ 

* 

★ 

★ 

it 

★ 


WELCOME  TO  THE  GERESS  ARRAY  DATA  CENTER  AT 
Ruhr-University  Bochiam 

ON-LINE  PROCESSING  DISPLAY  MANAGER 
(Version  3.4,  8.  Jun  94;  Dr.  M.  L.  Jost) 


it 

* 

★ 

★ 

■k 

★ 


***★**★★★★*★**★*★★**★★★★**★★****★*★★★★***★***★**★******♦*************** 


+  Lists  of  the  GERESS  on-line  processing  (extended  epi)  can  be  + 
+  sent  daily  to  you.  Please  give  us  your  e-mail  address  + 


TODAYS  DATE:  94159 

GERESS  **************************************************************** 


Cmds :  det 

GERESS  detections 

deti 

info  on  GERESS 

detector 

fk 

GERESS  f-k  results 

down 

info  on  GERESS 

downtimes 

epi 

GERESS  unrev.  locations 

bull 

GERESS  automat. 

.  Bulletin 

loc 

GERESS  rev.  locations 

tel 

GERESS  rev.  teles,  onsets 

plot 

TEKTRONIX  plot  module 

ALERTS  **************************************************************** 

nau 

NEIC  auto,  locations 

fei 

FEIS  alerts 

nar 

NEIC  rew.  auto,  locations 

sal 

SED  alerts 

ale 

NEIC  alerts 

ing 

ING  alerts 

ems 

EMSC  alerts 

nor 

NORSAR  alerts 

ONSETS  **************************************************************** 

grn 

GRSN  onsets 

yka 

YKA  onsets 

aus 

Austrian  Net  onsets 

vra 

VRAC  onsets 

csd 

KHC  PRU  onsets 

BULLETINS 

************************************************************* 

qed 

QEDs  from  NEIC 

Idr 

LDG  reg.  bull. 

(France) 

pde 

weekly  EDRs  from  NEIC 

Idt 

LDG  tel.  bull. 

(France) 

edr 

monthly  EDRs  from  NEIC 

sed 

SED  bulletin  (Switzerl) 

idc 

GSETT-3  EIDC  Bulletin 

igr 

IGN  reg.  bull. 

(Spain) 

ims 

IMS  Bulletin 

igt 

IGN  tel.  bull. 

(Spain) 

gol 

German  events  823-1989 

sol 

Swiss  events  80 

-87 

mur 

Mururoa  events 

isr 

Israel  events  07-93 

MOMENT  TENSORS  ******************************************************** 

cmt  CMTs  from  Harvard  sip  MTs  from  OSGS  (Sipkin) 

eri  CMTs  from  Tokyo-Univ. 


SERVICE  *************************************************************** 

on  lASPEI  onset  times  for  YODR  station 

sr  Seismotectonic  region  in  Central  Europe  for  given  event 

fie  Flinn-Engdahl  region  for  given  event 

doy  list  of  days  in  year 

help,  h  list  of  commands 

info  short  info  on  display  manager 

mess  short  message  to  GERESS/DISPLAY  manager 

exit,  X  logout,  end  session 

**************************************************v3.4**  8. Jun. 94# jost# 
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2.  MONITORING  NUCLEAR  TEST  SITES 


WITH  GERESS 


M.  L.  lost,  J.  Schweitzer,  and  H.-P.  Harjes 


2.0  ABSTRACT 


The  GERESS  array  in  south-easton  Gennany  has  been  proposed  to  parti¬ 
cipate  as  an  a-station  during  the  upcotning  GSETT-3  test  With  the  aim  of 
monitoring  a  comprehensive  test  ban  treaty,  special  interest  has  emerged  to 
asses  the  capability  of  this  new  array  to  monitor  the  known  nuclear  test  sites. 
During  the  GERESS  site  survey  (1988-1989)  and  4  1/2  years  of  GERESS 
operation  (Jan.  1990  -  June  1994),  44  presumed  nuclear  events  have  been 
recorded;  18  from  Nevada,  14  from  the  Tuamotu  Archipelago,  6  from  eastern 
Kazakhstan,  5  from  Lop  Nor,  and  1  from  Novaya  Zcmlya.  The  detection  thres¬ 
holds  at  GERESS  are  m,,  =  4.3  for  Nevada,  m^  =  4,0  for  Tuamotu  Archipelago, 
mb  =  4.6  for  Lop  Nor,  mb  =  4.3  for  Novaya  Zemlya,  and  nOb  =  4.5  for  eastern 
Kazakhstan.  For  Nevada  and  the  Tuamotu  Archipelago,  tiiese  detection  thres¬ 
holds  at  GERESS  correspond  to  a  yield  estimate  in  the  order  of  1  kt  at  each 
site.  These  path  specific  detection  thresholds  deviate  substantially  from  those 
derived  from  global  amplitude-distance  relations  (i.e.,  Gutenberg-Richter). 
Based  on  a  cross-correlation  method  and  a  cluster  analysis,  GERESS  records 
can  be  used  to  identify  specific  test  sites. 


2.1  INTRODUCTION 


The  Group  of  Scientific  Experts  (GSE)  of  the  Conference  of  Disarmament 
(CD)  in  Geneva  has  started  developing  a  final  concept  for  an  international 
seismic  data  exchange  system  (CD/1211,  1993)  providing  technical  support  for 
monitoring  a  comprehensive  nuclear  test  ban  treaty  (CTBT).  The  near  real 
time  system  includes  a  three-tiered  global  network  of  seismic  stations  where  the 
first  tier,  the  alpha-network,  consists  about  half  of  arrays  and  half  of  three- 
component  stations.  These  alpha-stations  are  designed  to  provide  not  only  the 
required  detection  threshold  tWughout  the  world  but  also  to  achieve  a  prelim¬ 
inary  event  location  which  is  then  improved  by  the  supplementary  data  of  the 
beta  and  gamma  networks.  In  this  respect,  the  a-network  of  about  60  globally 
distributed  stations  is  regarded  as  a  teleseismic  monitoring  system. 
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Recently,  the  GERman  Experimental  Seismic  System  (GERESS)  has  been 
proposed  to  participate  m  the  upcoming  GSE  TechMcal  Test  3  (GSETT-3)  as 
an  a-station.  Hence,  the  ct^abilities  of  the  GERESS  array  to  monitor  the 
presently  known  nuclear  test  sites  has  become  of  interest. 

Based  on  4  1/2  years  of  GERESS  data  (Jan.  1990  -  June  1994),  this  paper 
addresses  the  detection  capabilities  for  presumed  nuclear  events.  We  also  added 
10  events  from  the  GERESS  site  surveying  phase  (1988-1989)  to  include  data 
from  the  joint  verification  experiments  (JVE)  and  fiom  the  Kazakhstan  test  site. 
Due  to  the  amount  of  data,  we  focus  on  events  at  the  Nevada  Test  site  (NTS) 
and  events  in  the  Tuamotu  Archipelago. 


2.2  THE  GERESS  ARRAY 

Since  the  establishment  of  the  Giafenberg  array  (GRF)  in  1975,  consider¬ 
able  efforts  have  been  made  in  Germany  to  use  array  data  for  monitoring 
underground  nuclear  explosions  (Haijes  and  Hanka,  1986).  After  it  became 
clear  that  high  frequency  signals  play  an  in^xntant  role  for  detection  and 
identification  of  small  underground  nuclear  tests  at  regional  distances  (e.g., 
Mykkeltveit  et  al.,  1990;  Ringdal,  1990),  a  25-element  short  period,  small  aper¬ 
ture  array  in  south-eastern  Germany  was  plaimed.  This  array  complements  the 
broad-band,  large  aperture  Grafenberg  array  which  had  not  been  designed  for 
detection  purposes.  The  design  goals  of  GERESS  included  the  lowering  of  the 
global  detection  threshold  in  addition  m  mcHUtoring  local  and  regional  seismi¬ 
city  in  Central  Europe  and  the  Mediterranean  earthquake  zone. 

The  siting  survey  for  this  new  array  began  in  1987  (Haijes,  1989).  During 
this  survey,  single  stations  were  set  up  near  the  localities  Gross-Lichtenberg 
(GRLB  48.838  N,  13.720  E),  Haidmuehle  (HAID  48.800  N,  13.786  E),  and 
Sulzberg  (SULZ)  at  the  present  GERESS  C2  location.  These  single  stations 
operated  event  triggered  from  April-September  1988.  From  December  1988  to 
April  1989,  a  9-element  test-array  was  in  operation  at  the  present  GERESS  site 
(Haijes,  1990a).  The  25-element  GERESS  array  was  installed  in  cooperation 
with  Southern  Methodist  University  (Dallas,  Texas)  and  Ruhr-University 
Bochum  and  records  are  available  fiom  Apr.  1990  (Haijes,  1990b). 

GERESS  is  located  in  the  Bohemian  Forest  in  south-eastern  Germany  near 
the  border  to  Austria  and  the  Czech  Republic  (Fig.  2-1).  The  array  consists  of 
25  stations  with  vertical-component  short-period  Teled3me  Geotech  GS-13 
instruments  sampled  at  40  Hz.  In  addition,  four  of  the  sites  include  horizontal 
component  instruments.  At  the  key  station  of  the  array,  GEC2  (48.84511  N, 
13.70156  E,  1132  m),  a  STS-2  broad-band  instrument  is  sampled  at  80  Hz. 
The  geometry  of  the  array  is  based  on  concentric  rings  providing  an  overall 
aperture  of  about  4  km  (Fig.  2-2).  Data  from  the  array  are  continuously 
transmitted  to  the  NORSAR  data  center  in  KjeUer  and  to  the  Institute  of  Geo¬ 
physics,  Ruhr-University  Bochum. 
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The  GERESS  data  center  at  Bochum  is  operating  an  experimental  on-line 
processing  system  (Jost  1992b),  which  uses  software  developed  at  NORSAR 
(Mykkeltveit  and  Bungum,  1984;  Fyen,  1987)  and  modified  in  Bochum.  The 
automatic  system  detects  teleseismic  phases  in  addition  to  locating  regional 
events  (Haijes  et  al.,  1993;  Jost,  1993).  Recently,  Haijes  et  al.  (1994) 
evaluated  the  contribution  of  GERESS  to  global  seismicity  using  the  NEIC  bul¬ 
letin  as  reference  and  concluded  that  GERESS  has  exceUent  teleseismic  detec¬ 
tion  capabilities. 

GERESS  is  located  in  the  teleseismic  distance  range  with  respect  to  all 
nuclear  test  sites  that  were  active  1988  -  June  1994  (Fig.  2-3).  Table  2-1  indi¬ 
cates  the  mean  locations  of  these  test  sites,  the  epicentral  distances  to  GERESS, 
and  the  back-azimuths.  This  distance  range  may  actually  be  no  drawback  for 
the  task  of  monitoring  a  low-yield  test  ban  treaty.  For  example,  GERESS 
automatically  detected  a  P-wave  signal  from  a  chemical  calibration  experiment 
of  1  kt  charge  at  NTS  (sometimes  called  non-proliferation  test  (NPT))  whereas 
Lg  waves  were  barely  visible  at  the  low  noise  station  LTX  at  regional  distance 
(13®,  Lawson,  pers.  comm.). 


2.3  NEVADA  TEST  SITE 


From  the  total  of  24  underground  nuclear  explosions  conducted  during  the 
investigated  time  period,  18  waveforms  are  available  at  GERESS  in  axldition  to 
the  aforementioned  chemical  NPT.  4  of  the  24  events  (1990/9/20  9:15  local, 
1990/9/27  11:02  local,  1992/4/30  9:30  local,  and  1992/6/19  VICTORIA) 
were  too  small  to  produce  any  observable  signal  at  stations  in  Central  Europe. 
During  2  tests  (METROPOLIS  and  DISTANT  ZENITH),  GERESS  was  out 
Table  2-2  shows  parameters  of  events  at  NTS  from  1988-1993.  The  two  1988 
events  COMSTOCK  and  KEARSARGE  (JVE)  were  recorded  during  the  site 
survey  of  GERESS  at  station  GRLB. 

Figure  2-4  shows  the  waveftmns  of  aU  observed  events  from  Nevada, 
sorted  according  to  magnitude  detomined  at  GERESS  (Table  2-2).  All  traces 
start  10  seconds  before  the  theoretical  onset  time  calculated  ftom  the  published 
Department  of  Energy  (DOE)  parameters  as  given  in  NEIC  using  the  IASP91 
(Kennett  and  Engdahl,  1991)  travel  time  tables.  For  the  2  events  in  1988,  sin¬ 
gle  traces  are  shown.  For  the  events  in  1990-1993,  the  optimum  GERESS  beam 
is  displayed.  In  Figures  2-4  -  2-6,  we  filtered  all  available  vertical  traces  with  a 
third  ordCT  butterworth  band-pass  filtCT  from  0.5  Hz  -  2.5  Hz,  0.7  Hz  -  2.0  Hz, 
and  1.0  Hz  -  3.0  Hz,  respectively,  before  forming  the  displayed  optimum  beam. 
Only  forward  filtering  is  used  throughout  this  paper.  The  parameters  (back- 
azimuth  and  slowness)  to  calculate  the  optimum  GERESS  beam  were  obtained 
from  fk  analysis  using  all  available  vertical  traces  at  GERESS. 

In  Figures  2-4  -  2-6,  we  clearly  see  a  scatter  of  +/-  1  sec  for  the  observed 
onset  times  at  GERESS.  A  first  explanation  would  be  that  the  absolute  timing 
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of  the  GERESS  array  was  not  stable  and  produced  these  deviations.  But  this 
conclusion  could  be  rejected  because  all  GERESS  onset  times  show  a  constant 
travel  time  difference  of  7.5  -  7.8  sec  with  respect  to  published  onset  times  of 
these  events  at  the  Grafenberg  array.  It  is  obvious  that  the  same  absolute  time 
scatter  could  not  occur  at  two  different  arrays.  In  addition,  we  verified  the  rela¬ 
tive  onset  times  of  MONTELLO  and  HOUSTON  at  stations  Berggiesshilbel 
(BRG),  Collmberg  (CLL),  and  Kaspersk’e  Hory  (KHC).  The  observed  scatter 
of  onset  times  can  be  caused  in  part  by  the  velocity  structure  at  NTS 
(McLaughlin  et  al.,  1987;  Ferguson  et  al.,  1989).  Because  a  time  dependent 
travel  time  between  NTS  and  the  two  arrays  in  Germany  is  not  realistic  our 
preferred  hypothesis  is  that  the  published  source  times  are  not  correct 

For  the  events  in  1990-1993,  we  measured  arrival-times,  amplitudes,  signal 
to  noise  ratios  (SNR),  and  periods  on  the  beam  traces  (Table  2-2).  It  is  not 
surprising,  that  very  small  slowness  values  (the  theoretical  value  firom  IASP91 
is  5.13  sec/deg)  were  determined,  since  the  GERESS  aperture  of  about  4  km 
limits  the  resolution  of  slowness  measurements  in  the  teleseismic  distance 
range. 

The  nuclear  tests  AUSTIN,  FLOYDATA,  GALENA,  and  DIVIDER  were 
not  detected  by  the  automatic  STA/LTA  detector  at  GERESS  but  show  weak 
signals  at  their  expected  onset  times  in  Figures  2-4  -  2-6.  The  measured  param¬ 
eters  (Table  2-2)  indicate  an  upper  limit  for  these  events.  In  the  filter-bmd  0.7 
-  2.0  Hz  (Fig.  2-5),  weak  onsets  of  AUSTIN,  GALENA,  and  DIVIDER  seem 
to  be  present.  Parameters  were  measured  on  these  beams  (theoretical  beam 
parameters).  Using  the  optimum  detection  filter-band  for  GERESS  (lost 
1992a),  i.e.,  1.0  -  3.0  Hz,  a  weak  onset  for  FLOYDATA  is  visible  (Fig.  2-6). 
Parameters  of  FLOYDATA  were  measured  on  the  corresponding  theoretical 
beam.  For  all  4  events,  the  fk  analyses  produced  less  reliable  results  due  to  the 
low  SNR. 

A  detection  threshold  of  mb  (NEIQ  =  4.3  is  estimated  using  16  NTS- 
events  recorded  at  GERESS  in  1990-1993.  This  is  done  by  comparing  the 
measiued  SNR  values  (Table  2-2)  with  the  threshold  of  the  GERESS  STA/LTA 
detector  of  3.8  on  a  test-beam  (e.  g.,  indirect  estimation  method;  Ringdal, 
1975).  The  calculated  GERESS  nib-values  are  in  good  agreement  witii  the  net¬ 
work  magnitudes  published  by  NEIC.  GERESS  did  not  detect  COSO  with 
nib  =  4.2  (GERESS).  Also  GALENA  nib  =  4.0  (GERESS)  was  only  detected 
after  reprocessing  at  a  SNR  of  2.1.  On  the  other  hand,  HUNTERS  TROPHY 
with  nib  =  4.2  (GERESS)  and  also  the  NPT  with  mb  =  4.0  (GERESS)  were 
clearly  detected  with  a  SNR  above  5.  The  influence  of  noise  (day  /  night;  wind) 
as  well  as  station  outages  affect  the  SNR  and  can  explain  this  variability. 

In  order  to  obtain  a  rough  estimate  on  the  yield  from  the  single  site 
GERESS,  we  fit  a  straight  line  through  the  2  known  yields:  The  NPT 
(nib  (GERESS)  =  4.0)  had  an  equivalent  yield  of  1  kt  following  Denny  and 
Zucca  (1993).  The  JVE  (KEARSARGE)  had  a  yield  close  to  150  kt  but  pro¬ 
duced  a  body  wave  magnitude  of  only  5.4  (GERESS).  On  the  other  hmd, 
explosions  like  JORNADA  produced  a  magnitude  of  5.9  at  a  comparable  yield 
(Gupta  et  al.,  1989).  JORNADA  was  detonated  below  the  water  table  and 
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KEARSARGE  above,  causing  less  coupling  for  the  latter.  Gupta  et  al.  (1989) 
found  that  teleseismic  mt,  values  for  shots  below  the  water  table  are  0.3  -  0.5 
magnitude  units  larger  than  for  shots  above  the  water  table  at  a  fixed  yield. 
Following  these  authors,  we  assumed  that  the  water  table  influence  is  about  0.4 
magnirnde  units,  scaling  our  observed  m|,  =  5.4  (GERESS)  for  KEARSARGE 
to  5.8  for  a  regular  150  kt  NTS  shot  below  the  water  table.  We  obtained 

logio  (yield)  =  1.21  mb  (GERESS)  -  4.84  .  (1) 

The  resulting  yield  estimates  from  GERESS  are  given  in  Table  2-3.  Fbr  com¬ 
parison,  magnitudes  and  yields  determined  at  GRF  by  Schlittenhardt  (1988l>- 
1992)  are  included.  Using  Nuttli’s  (1986,  1987)  linear  relation  between  1% 
(Lg)  and  yield  and  taking  proper  care  of  the  magnitude  bias  between  P-wave 
and  Lg-wave  magnitudes,  we  obtained  further  yield  estimates  fit>m  GERESS 
P-wave  magnitudes  (Table  2-3).  Recently,  Goldstein  et  al.  (1994)  estimated  the 
yield  for  BRISTOL  as  10  kt.  They  argued  that  their  value  could  be  30  %  lower 
due  to  differences  in  coupling.  We  estimated  for  BRISTOL  7  kt  from  (1),  and 
8  kt  using  Nuttli’s  relation.  Essentially,  the  yield  estimates  from  (1),  Nuttli’s 
relation,  and  GRF  are  not  significantly  different  considering  the  lack  of  detailed 
shot  information  (e.g.,  water  table,  spall).  In  addition,  effects  of  local  geology 
can  affect  teleseismic  P-wave  magnitudes  introducing  uncertainQr  to  the  yield 
estimate  (McLaughlin  et  aL,  1987;  Cormier,  1987).  In  Table  2-3,  a  yield  esti¬ 
mation  was  calculated  for  all  events.  This  estimation  together  with  the 
automatic  detection  of  the  NPT  suggest  that  GERESS  can  detect  events  in  the 
order  of  1  kt  at  NTS. 

Furthermore,  Table  2-3  suggests  that  GERESS  can  observe  nuclear  tests  at 
NTS  clearly  below  a  yield  estimate  of  1  kL  Although  small  events  like  FLOY- 
DATA  or  GALENA  were  not  automatically  detected,  their  visible  waveforms 
can  be  used  to  confirm  a  location  hypothesis  within  the  firamework  of  GSETT- 
3. 

In  order  to  identify  subgroups  of  the  Nevada  data  set,  we  performed 
cross-correlations  of  the  optimum  beams.  Figure  2-7  shows  the  correlation 
matrix  obtained  by  using  Ae  maximum  cross-correlation  value  for  the  filter- 
band  0.5-2.5  Hz  and  the  event  numbers  from  Table  2-3.  We  clearly  see  that 
the  small  events  (mb  <  4.3;  events  #  1  -  #  7)  do  not  have  a  correlation  larger 
than  0.54,  possibly  due  to  the  small  SNR.  On  the  other  hand,  the  cross- 
correlation  sequences  for  FLOYDATA,  AUSTIN,  DIVIDER,  and  GALENA 
indicate  that  Ae  observed  waveforms  are  from  the  corresponding  events  at 
NTS.  For  the  larger  events  (mb  >  4.6;  events  #  8  -  #  19),  the  correlation 
reaches  values  of  up  to  0.91. 

We  further  performed  a  single  linkage  clustering  analysis  (Everitt,  1993). 
Figure  2-8  shows  the  corresponding  denthogram.  We  see  that  JUNCTION  (# 
17),  MONTELLO  (#  10),  HOUSTON  (#  14),  and  COMSTOCK  (#  13)  form  a 
cluster.  Also  KEARSARGE  (#  12)  and  BEXAR  (#  15)  form  a  subgroup.  An 
explanation  for  the  latter  could  be  the  nearly  identical  shot-locations  (Fig.  2-9). 
Also  HUNTERS  TROPHY  (#  6)  and  the  NPT  (#  5)  appear  somewhat  related, 
although  the  significance  level  is  rather  low.  Obviously,  the  shots  at  Pahute 
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Mesa  and  Yucca  Flat  do  not  form  independent  clusters.  This  result  agrees  with 
previous  observations  by  Murphy  et  al.  (1989)  who  found  that  the  seismic 
source  characteristics  of  both  testing  areas  are  identical  for  tests  below  the 
water  table.  In  general,  we  observe  large  correlation  coefficients  for  events  at 
Pahute  Mesa  (Fig.  2-9).  There  the  depth  of  burial  is  larger  compared  to  events 
at  Yucca  Flat  inplying  a  larger  yield.  The  lower  correlation  values  for  most  of 
the  shots  at  Yucca  Rat  are  probably  due  to  their  small  size  and  hence  the  small 
SNR  of  the  beams  entering  the  correlation  analysis.  Another  reason  may  be  the 
variability  of  the  geologic  structure  at  Yucca  Rat  (e.g.,  McLaughlin  et  al., 
1987;  Cormier,  1987;  Ferguson,  1988). 


2.4  TUAMOTU  ARCHIPELAGO 


In  the  epicentral  distance  range  between  143°  and  147°,  the  structure  of 
the  Earth’s  cote  amplifies  -  like  a  focusing  lens  -  seismic  waves  which  travel 
through  the  Earth’s  core  as  PKP-phases.  The  comparison  of  the  daily  detection 
lists  of  GERESS  with  bulletins  of  international  seismic  data  centers  shows  that 
GERESS  is  therefore  very  sensitive  for  the  seismicity  in  the  southern  Pacific. 
The  two  French  nuclear  test  sites  Mururoa  and  Fangataufa  are  located  very 
close  to  each  other  in  the  Tuamotu  Archipelago  (Table  2-1).  They  have  an  epi¬ 
central  distance  to  GERESS  of  about  145.5°,  right  where  the  maximum  of  the 
PKP-caustic  is  observed.  Therefore,  we  expect  a  low  detection  level  for  seismic 
events  at  these  test  sites. 

The  waveforms  of  14  presumed  French  nuclear  tests  in  1988-1991  are 
shown  in  Figure  2-10  (optimum  beams,  single  traces  for  events  in  1988,  all 
filtered  0.5  -  2.5  Hz).  The  events  are  ordered  with  respect  to  increasing  yield 
estimate  (Table  2-4).  Table  2-4  gives  the  corresponding  parameters  measured 
on  the  optimum  beams  (events  in  1990  and  1991)  or  on  the  single  traces  for  the 
events  in  1988  (recorded  at  station  HAID).  The  French  authorities  have  not  yet 
published  any  parameters  of  their  tests.  Therefore,  the  start  time  of  the 
waveform  for  the  event  on  Nov.  14,  1990  was  chosen  to  10  seconds  before  the 
measured  onset  time.  The  other  traces  woe  aligned  to  its  first  minimum.  In 
1991,  one  of  the  presumed  explosions  was  not  detected  by  GERESS  (Jul  5, 
1991;  18:00).  This  event  was  very  small  and  was  only  detected  by  a  T-phase 
observation  at  Raratonga  (RAR).  The  yield  was  estimated  by  DSIR  (Smith, 
pers.  communication)  as  less  than  0.3  kt  Rgure  2-10  also  displays  the 
GERESS  theoretical  beam  at  the  expected  onset  time  for  this  event,  but  no  con¬ 
clusive  signal  could  be  identified. 

The  body  wave  magnitude  mi,  is  not  defined  for  events  with  an  epicentral 
distance  of  more  than  about  100°  and  no  direa  m|,  determination  is  possible  for 
these  events  from  GERESS.  Schlittenhardt  (1988a)  introduced  a  body  wave 
magnitude  for  PKP-observations  of  the  French  explosions  at  the  Grafenberg 
array.  Although  the  distance  between  the  reference  station  Al  of  GRF  and 
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GERESS  is  only  about  200  Imi,  we  cannot  use  his  results  to  estimate  m{,-values 
because  the  site  conditions  differ  substantially.  The  amplitude-distance 
behavior  of  the  PKP-phases  varies  drastically  in  the  distance  range  near  145.5° 
(e.  g.,  Houaid  et  al.,  1993),  and  we  COTrected  the  observed  log  (A/T)  values 
with  the  attenuation  curve  of  Blandford  and  Sweetser  (1973)  for  PKP-phases. 
With  the  correction  value  of  3.50  for  an  epicentral  distance  of  145.5°,  we  deter¬ 
mined  the  mb(PKP)-values  for  the  presumed  French  tests.  For  11  of  these 
events,  our  mj,(PKP)-values  can  be  compared  with  mi,-determinations  by 
NEIC.  Althougjii  our  mb  values  are  slightly  smaller  than  the  netwoik  magni¬ 
tudes  by  NEIC,  the  general  agreement  between  the  nib-values  calculated  with 
different  methods  confirms  the  proposed  usage  of  PKP-phases  for  regular 
world-wide  monitoring  (Haijes,  1985). 

Figure  2-11  shows  the  estimated  mb(PKP)-values  for  GERESS  plotted 
against  the  yields  estimated  by  DSIR  (Table  2-5,  Smith,  pers.  communication). 
The  line  follows  the  linear  relationship  (least  squares  fit)  between  log(yield)  and 
mb  (PKP)-values; 

logic  (yield)  =  13i  nOb  (PKP)  -  4.88  .  (2) 

This  relation  is  very  similar  to  the  results  of  Schlittenhardt  (1988a)  who  got  for 
mb(PKP)  based  on  GRF  measurements: 

logic  (yield)  =  1.12  mb  (PKP)grf  -  4.58  .  (3) 

The  two  solid  lines  in  Figure  2-11  give  the  upper  limits  for  the  smallest  not 
observed  explosion.  On  average,  the  GERESS  detection  threshold  is  about 
mb(PKP)  =  3.8  corresponding  to  a  yield  estimate  of  0.7  kL  For  the  network 
(NEIC)  Tnagnitiidft,  the  GERESS  detection  threshold  is  4.0. 

Nuclear  tests  in  the  Tuamotu  Archipelago  have  been  repented  to  be 
detonated  either  on  the  island  of  Mururoa  or  on  the  island  of  Fangataufa  with 
about  40  km  in  between.  The  corresponding  variation  in  slowness  and  azimuth 
for  PKP-onsets  is  so  small  that  even  large  aperture  arrays  (e.g.,  GRF)  are  not 
able  to  discriminate  between  the  2  test  sites.  Therefore,  we  used  a  different 
approach  and  performed  a  cross-correlation  of  all  recorded  events  (1988-1991). 
We  used  a  30  second  time  window  of  the  optimum  GERESS  beams  filtered 
between  0.5-2.5  Hz.  Figure  2-12  shows  a  cross-correlation  sequence  of  the 
event  on  Nov.  21,  1990  with  all  events  in  1990-1991.  Note  that  the  cross¬ 
correlation  coefficients  of  the  three  events  on  the  bottom  of  the  figure,  report¬ 
edly  placed  on  Fangataufa,  are  significantly  different  (very  small)  firom  those 
reportedly  placed  on  Mururoa.  The  cross-correlation  coefficients  are  also  small 
for  the  first  two  events.  The  first  event  has  been  classified  as  not  observed,  and 
the  second  (yield  in  the  order  of  1  kt),  although  observed,  is  contaminated  with 
too  much  noise.  On  the  other  hand,  the  cross-correlation  sequence  indicates 
that  the  observed  waveform  is  from  the  corresponding  event  at  Tuamotu.  Fig¬ 
ure  2-13  shows  the  cross-correlation  matrix  obtained  by  using  the  maximum 
cross-correlation  value  in  the  filter-band  0.5-2.5  Hz  and  the  event  numbers  from 
Table  2-5.  Table  2-6  lists  the  minimum  and  maximum  correlation  coefficients 
for  each  event  with  respect  to  the  Mururoa  and  Fangataufa  clusters.  Repeating 
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the  correlation  analysis  including  all  recorded  events  in  1988  shows  that  the 
difference  of  correlation  coefficients  between  those  two  clusters  decreases.  We 
explain  this  by  the  larger  noise  level  of  the  single  station  records  (events  in 
1988).  Other  pass-bands  (0.7-2.0  Hz,  1.0-3.0  Hz)  generally  support  the  present 
finding,  although  the  observed  grouping  of  correlation-coefficients  is  optimum 
for  most  events  in  the  pass-band  of  0.S-2.S  Hz.  The  cross-ccnrelation  sequence 
of  the  event  on  May  25,  1988  (#  11)  which  had  the  largest  reported  yield  on 
Mururoa  (83  kt)  also  agrees  with  the  present  finding,  where  an  unexplained  de- 
correlation  with  the  event  on  June  2,  1990  (#  8)  is  observed  (Fig.  2-13).  Figure 
2-14  shows  the  dendrogram  based  on  a  single  linkage  clustering  analysis.  We 
clearly  see  the  different  Fangataufa  and  Mururoa  clusters.  The  event  #  2 
(07/05/1991;  1  kt)  has  a  very  small  SNR  which  makes  a  reasonable  analysis 
impossible.  On  the  other  hand,  we  see  some  variability  in  the  Mururoa  cluster, 
e.g.  events  #  4  (18/5/1991),  #  10  (21/11/1990),  and  #  9  (15/07/1991)  appear  to 
form  a  subgroup  and  also  the  events  #  6  (11/05/1988),  #  5  (04/07/1990),  and  # 
7  (14/06/1991).  This  observation  may  result  from  the  observed  local  site 
effects  (Crusem  and  C!aristan,  1992). 

Following  Guille  et  al.  (1993),  the  atoll  of  Mururoa  (28  km  by  10  km;  age 
11.8-10.6  My)  is  part  of  the  volcanic  chain  from  the  hot-spot  near  the  island  of 
Pitcairn.  Mururoa  sits  on  top  of  a  submarine  plateau  of  130  km  by  30  km 
oriented  N80E.  The  existence  of  this  plateau  is  intimately  linked  to  the  Austral 
fracture  zone  which  explains  the  atypical  mientation  of  this  plateau.  Although 
the  island  of  Fangataufa  (10  km  by  5  km;  age  11.5-9.6  My)  is  situated  only  40 
km  south-east  of  Mururoa,  the  influence  of  this  fracture  Zone  is  minor.  The 
morphology  of  the  basement  is  star  shaped,  quite  different  from  the  elongated 
plateau  of  Mururoa.  Typical  tectonic  directions  on  Fangataufa  follow  those  of 
the  southern  Pacific  (NS,  N80E,  and  N115E).  The  bathymetry  also  supports  a 
generic  difference  between  both  islands:  the  strait  separating  both  islands  shows 
water-depths  larger  than  3  km  fen-  about  20  km.  As  a  result,  both  islands  exhi¬ 
bit  different  tectonic  settings  that  can  be  responsible  for  the  observed  clustering. 

We  conclude  that  based  on  the  limited  data  set  (14  events),  the  optimum 
GERESS  beams  can  be  used  in  a  correlation  analysis  (waveform  matching)  to 
successfully  separate  the  Mururoa  and  Fangataufa  clusters.  French  authmities 
have  not  yet  published  any  information  on  their  tests.  Therefore  we  cannot 
know  to  which  degree  the  above  results  reflect  different  device  types. 
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2.5  LOP  NOR 


GERESS  lecoided  5  of  the  6  jnesumed  nuclear  explosions  at  the  Chinese 
test  site  Lop  Nor  in  1990-1994.  Figure  2-15  displays  the  optimum  beams 
filtered  0.5  -  2.5  Hz.  Parameters  measured  on  the  optimum  beams  are  given  in 
Table  2-7.  The  observed  body  wave  magnitudes  are  unexpected  low  compared 
to  the  network  magnitudes  (NEIC).  This  observed  discrepancy  suggests  that 
GERESS  has  a  significant  station  residual  for  Nop  Lor  events.  However, 
GERESS  generally  shows  a  negative  magnitude  residual  for  events  in  western 
China  and  the  Ifimalaya  region  of  -0.51  +/-  0.26  magnitude  units.  With  a 
STA/LTA  detection  threshold  of  3.8  of  the  GERESS  detector,  this  implies  a 
detection  threshold  for  events  at  Lop  Nor  of  mj,  (NEIC)  =  4.6. 

For  4  events,  surface  wave  magnitudes  could  be  measured  at  GERESS 
(Table  2-7).  However  for  the  events  in  1990  and  1993,  the  magnitude  values 
give  an  upper  bound  only  since  the  traces  were  disturbed  by  interfering  noise. 
The  classical  mi,  /  MS-ratios  at  GERESS  clearly  confirm  the  explosion  source 
type  for  all  4  events.  On  the  other  hand.  Lop  Nor  recmds  show  a  compara¬ 
tively  strong  Love-wave  radiation  towards  European  stations,  which  is  not 
expected  for  an  explosion  type  source.  Note  that  die  surface  wave  magnitude 
measured  on  the  horizontal  component  instruments  is  significantly  larger  than 
the  magnitude  measured  on  the  vertical  instrument  (Table  2-7)  because  the 
Love  wave  shows  the  largest  anoplitude  in  the  surface  wave  group. 

The  test  on  May  21,  1992  was  the  largest  underground  explosion  which 
has  ever  been  observed  with  broad-band  stations  world-wide.  To  demonstrate 
the  simple  pulse  form  of  relatively  high  frequency  for  this  large  event.  Figure 
2-16  shows  the  deconvolved  P  and  PcP  onsets,  fii  both  cases,  a  long  periodic 
down  swing  after  the  first  onset  is  seen  which  can  be  addressed  as  an  effect  of 
tectonic  release  related  with  this  event  (Gao  and  Wallace,  1992). 


2.6  NOVAYAZEMLYA 


On  October  24,  1990,  the  Soviet  Union  carried  out  its  last  nuclear  test 
(Table  2-8).  This  test  was  conducted  on  Novaya  Zemlya  at  a  distance  of  30.4° 
to  GERESS.  The  records  are  dominated  by  the  very  prominent  P-wave  coda 
(Fig.  2-17,  optimum  beam,  filtered  0.5  -  2.5  Hz),  and  a  very  complex  pulse  of 
the  direct  P-wave.  This  feature  must  be  explained  by  contributions  from  the 
direct  waves  of  the  upper  mantle  triplications,  by  converted  phases,  multipath- 
ing,  and  by  a  lateral  heterogeneous  structure  of  crust  and  upper  mantle  in 
Europe.  The  observed  magnitude  mi,  =  5.5  for  this  explosion  is  in  agreement 
with  NEIC  (mi,  =  5.7).  The  classical  mi,  /  MS-ratio  (5.5  /  4.4  =  1.3)  at  GERESS 
clearly  confirms  the  explosion  source.  With  the  observed  SNR  of  90.3  the 
detection  threshold  at  GERESS  can  be  estimated  at  mi,  (NEIC)  =  4.3. 
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Using  Nuttli’s  (1986,  1988)  linear  relation  between  m{,  (Lg)  and  yield  and 
taking  proper  care  of  the  magnitude  bias  between  P-wave  and  Lg-wave  magni¬ 
tudes,  we  obtained  a  yield  estimate  of  50  kt  from  the  GERESS  P-wave  magni¬ 
tude.  Schlittenhardt  (1988b-1992)  determined  the  yield  for  the  Novaya  Zemlya 
test  as  70  kt  based  on  GRF  records  (mb(GRF)=5.9).  For  GERESS,  we 
estimated  a  detection  threshold  in  the  order  of  3  kt  Note  that  this  value  is 
based  on  one  single  shot  only. 


2.7  EASTERN  KAZAKHSTAN 


During  the  site  survey  for  GERESS,  6  events  from  the  former  Soviet  test 
site  near  the  city  of  Semipalatinsk  were  recorded.  All  tests  were  conducted  at 
the  Shagan  River  (Balapan)  subsection  of  the  testing  area  (e.g.,  Thurber  et  al., 
1993).  We  recorded  the  event  on  May  4,  1988  at  the  temporary  station  HAID, 
and  the  tests  on  Jun  14  and  Sep  14  (JVE),  1988  at  the  temporary  station  GRLB 
(Haijes,  1989).  The  three  events  in  Dec.  1988,  Jan.  and  Feb.  1989  were 
recorded  with  the  temporary  9-element  GERESS  test  array  (Harjes  1990a). 
Figure  2-18  shows  the  single  traces  and  the  optimum  beams,  req>ectively, 
filtered  between  0.5  Hz  -  2.5  Hz.  Except  for  the  event  on  June  14,  1988,  mag¬ 
nitudes  measured  at  the  GERESS  site  correlate  well  with  the  ISC  netwmk  mag¬ 
nitudes  (Table  2-9).  Based  on  these  6  events,  an  average  detection  threshold  of 
m{,  (NEIC)  =  4.7  has  been  determined.  This  estimate  gives  an  upper  bound 
only,  since  the  25-element  GERESS  array  shows  a  significandy  better  SNR 
improvement  by  beamforming  than  a  single  station  or  the  9-element  test  array. 
Hence,  the  detection  threshold  for  the  25-element  GERESS  array  can  be 
estimated  to  be  lower  than  4.5  for  eastern  Kazakhstan.  This  estimate  is  conser¬ 
vative  because  the  comparison  of  the  technical  installation  of  GERESS  and  the 
temporary  test  array  suggests  a  sensitivity  surpassing  the  expected  Vn  SNR 
improvement  Ringdal  (1990)  determined  the  50%  detection  threshold  for 
NORESS  as  m^  =  2.7  -  3.3.  The  discrepancy  to  die  estimate  for  GERESS  is  in 
part  due  to  the  well  dociunented  focusing  effects  in  the  upper  mande  below  the 
test  site  towards  Scandinavia.  Barker  and  Murphy  (1992)  determined  station 
residuals  in  Europe  which  amount  to  more  than  1  magnitude  unit 

In  Table  2-10,  we  compare  various  yield  determinations  for  the  6 
presumed  tests.  Using  Nuttli’s  (1986,  1987)  linear  relation  between  rUb  (Lg) 
and  yield  and  taking  proper  care  of  the  magnitude  bias  between  P-wave  and 
Lg-wave  magnitudes,  we  obtained  yield  estimates  from  GERESS  P-wave  mag¬ 
nitudes.  Included  are  also  yield  estimates  by  Schlittenhardt  (1988b-1992)  from 
GRF  records.  Based  on  the  Lg  wave  at  NORSAR,  Hansen  et  al.  (1990)  calcu¬ 
lated  a  yield  of  110  kt  for  the  JVE  on  Sep.  14,  1988.  This  value  agrees  with 
our  estimate  of  108  kt  at  GERESS.  From  hydrodynamic  measurements,  the 
yield  was  determined  from  US  and  USSR  scientists  at  115  kt  and  122  kt, 
respectively  (Sykes  and  Ekstrom,  1989).  Ringdal  et  al.  (1992)  and  also 
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Ekstrom  and  Richards  (1994)  estimated  source  parameters.  Considering  the 
variability  in  published  yield  determinations,  our  single  station  value  can  give  a 
valid  first  order  estimate  on  yield.  Hence,  for  the  25-element  GERESS  array,  a 
detection  threshold  in  the  order  of  2  kt  is  estimated. 


2.8  CONCLUSIONS 

During  the  GERESS  site  survey  (1988-1989)  and  4  1/2  years  of  GERESS 
operation  (Jan.  1990  -  June  1994),  44  presumed  nuclear  tests  have  been 
recorded:  18  from  Nevada,  14  from  the  Tuamotu  Archipelago,  6  fix>m  eastern 
Kazakhstan,  5  from  Lop  Nor,  and  1  from  Novaya  Zemlya.  The  detection  thres¬ 
holds  at  GERESS  are  estimated  as  mi,  =  4.3  for  Nevada,  noi,  =  4.0  for  Tuamotu 
Archipelago,  lUi,  =  4.6  for  Lop  Nor,  mi,  =  4.3  for  Novaya  Zemlya,  and  mi,  = 
4.5  for  eastern  Kazakhstan.  These  path  specific  detection  thresholds  deviate 
substantially  from  those  derived  firom  global  amplitude-distance  relations  (i.e., 
Gutenberg-Richter).  They  can  also  not  be  explained  with  a  simple  station  bias 
because  GERESS  generally  shows  a  negative  station  residual  of  0.6  magnitude 
units  +/-  0.35  (Haijes  et  al.,  1994).  Based  on  2  known  yields,  a  magnitude 
yield  relation  was  calculated  for  Nevada,  indicating  that  events  with  a  yield 
estimate  in  the  order  of  1  kt  can  be  detected  at  GERESS.  Based  on  yield  esti¬ 
mates  by  DSIR  (Smith,  pers.  comm.),  a  magnitude  yield  relation  was  deter¬ 
mined  for  events  in  the  Tuamotu  Archipelago.  The  detection  threshold  of  mi,  = 
4.0  would  correspond  to  events  with  a  yield  estimate  in  the  order  of  0.7  kt 
However,  for  Novaya  Zemlya  and  the  Kazakhstan  test  sites,  yields  have  been 
determined  by  various  authors.  Although  the  data  quantity  and  quality  at 
GERESS  for  these  test  sites  is  not  ideal,  we  estimated  a  detection  threshold  of 
3  kt  for  Novaya  Zemlya  and  2  kt  for  Kazakhstan.  Based  on  a  cross-correlation 
method  and  a  cluster  analysis,  GERESS  records  can  be  used  to  discriminate 
between  the  Tuamotu  and  Fangataufa  events. 

This  case  study  also  shows  that  caution  should  be  exercised  by  using 
results  of  theoretical  studies  on  network  detection  c^abilities.  Instead,  path 
specific  station  detection  thresholds  should  be  used  for  a  fixed  station  network 
to  derive  reliable  network  detection  capabilities. 

A  reliable  prediction  of  the  performance  of  a  future  global  monitoring  net¬ 
work  would  greatly  benefit  firom  information  about  exact  parameters  of  previ¬ 
ous  nuclear  tests.  In  this  respect,  the  Chinese  and  French  authorities  should  be 
appealed  to  release  the  corresponding  information  even  without  exact  yield 
information. 
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Table  2-1:  Distances  of  Nuclear  Test  Sites  to  GERESS 


Country 

Test  Site 

Lat 

Long 

A 

[deg] 

BAZ 

[deg] 

China 

Lop  Nor 

41.6  N 

88.7  E 

51.3 

68 

France 

Fangataufa 

22.2  S 

138.8  W 

145.6 

311 

France 

Mururoa 

21.9  S 

139.0  W 

145.4 

311 

Great  Britain 

Nevada 

37.2  N 

116.2  W 

83.5 

322 

USA 

Nevada 

37.2  N 

116.2  W 

83.5 

322 

USSR 

Kazakhstan 

49.9  N 

78.8  E 

41.2 

63 

USSR 

Novaya  Zemlya 

73.4  N 

54.7  E 

30.4 
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Table  2-2:  Parameters  of  Events  from  Nevada  Test  Site  Recorded  at  GERESS  1988-1993 
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1:  Single  station  of  GERESS  site  investigation 
2:  GERESS  out 

3:  Onset  at  GERESS  very  weak,  parameters  indicate  upper  threshold 
4:  Amplitude  measured  on  trace  filtered  1. 0-3.0  Hz 


Table  2-3:  Magnitude  Yield  Relations  for  NTS 


Yr 

Mo 

Day 

Name  /  Site 

“b 

GRF 

Yield 

GRF 

“b 

GERESS 

Yield 

GERESS 

Yield 

Nuttli 

No. 

1988 

Jun 

2 

Comstock 

5.4 

60 

5.4 

49 

68 

13 

1988 

Aug 

17 

Kearsarge,  JVE 

5.5 

75 

5.4 

150 

68 

12 

1990 

Jun 

13 

Bullion 

5.8 

150 

5.7 

114 

150 

19 

1990 

Jun 

21 

Austin 

<3.9 

<0.8 

2 

1990 

Jul 

25 

Mineral  Quarry 

5.0 

20 

4.7 

7 

8 

9 

1990 

Oct 

12 

Tenabo 

5.6 

90 

5.5 

65 

92 

16 

1990 

Nov 

14 

Houston 

5.5 

70 

5.4 

49 

68 

14 

1991 

Mar 

08 

Coso 

4.2 

2 

7 

1991 

Apr 

04 

Bexar 

5.7 

110 

5.5 

65 

92 

15 

1991 

Apr 

16 

Montello 

5.4 

50 

5.3 

37 

50 

10 

1991 

Aug 

15 

Floydata 

<3.6 

<0.3 

1 

1991 

Sep 

14 

Hoya 

5.6 

90 

5.6 

86 

125 

18 

1991 

Oct 

18 

Lubbock 

5.2 

30 

5.4 

49 

68 

11 

1991 

Nov 

26 

Bristol 

4.9 

10 

4.7 

7 

8 

8 

1992 

Mar 

26 

Junction 

5.6 

100 

5.6 

86 

125 

17 

1992 

Jun 

23 

Galena 

<4.0 

<1 

4 

1992 

Sep 

18 

Hunters  Trophy 

4.2 

2 

6 

1992 

Sep 

23 

Divider 

<3.9 

<0.7 

3 

1993 

Sep 

22 

NPT 

4.0 

1 

5 

The  GERESS  yields  are  from  log  y  =  1.21  mb  -  4.84.  The  Nuttli  yields  were  calculated  using 
his  1986  paper  and  GERESS  magnitudes  (see  text). 

The  No.  column  indicates  a  sorting  of  the  events  according  to  mb(GERESS). 
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Table  2-4:  Parameters  of  Events  from  Tuamotu  Archipelago  Recorded  at  GERESS  1988-1991 
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F:  event  probably  on  Fangataufa 

1:  Single  station  of  GERESS  site  investigation 

2:  event  not  observed,  parameters  measured  on  noise  at  theor.  onset  time 


Table  2-5:  Magnitude  Yield  Relations  for  Tuamotu  Archipelago 


Yr 

Mo 

Day 

Name  /  Site 

mb 

GRF 

Yield 

GRF 

mb 

GERESS 

Yield 

DSIR 

No 

1988 

05 

11 

Muniroa 

5.3 

6 

1988 

05 

25 

Mururoa 

65 

5.3 

83 

11 

1990 

06 

02 

Mururoa 

5.4 

40 

5.1 

8 

1990 

06 

07 

Mururoa 

4.8 

4.3 

3 

3 

1990 

06 

26 

Fangataufa 

5.7 

5.5 

12 

1990 

07 

04 

Mururoa 

5.2 

4.7 

18 

5 

1990 

11 

14 

Fangataufa 

5.7 

120 

5.5 

117 

14 

1990 

21 

Mururoa 

5.7 

35 

5.4 

36 

1991 

05 

07 

Mururoa 

4.6 

4 

1 

2 

1991 

05 

18 

Mururoa 

5.4 

16 

4 

1991 

05 

29 

Fangataufa 

5.8 

5.4 

13 

1991 

06 

14 

Mururoa 

5.5 

5.1 

28 

7 

1991 

07 

05 

Mururoa 

<0.5 

<3.3 

<0.3 

1 

1991 

07 

15 

Mururoa 

5.6 

70 

5.2 

34 

9 

The  No.  column  indicates  a  sorting  of  the  events  according  to  yield  estimate  (DSIR). 
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Table  2-6:  Results  of  Correlation  Analysis  of  Events  from  Tuamotu  Archipelago 
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(1) :  This  event  decorrelates  with  that  one  on  June  2,  1990  (cc=0.40) 

(2) :  Difference  in  correlation  coefficients  is  0.13  in  pass-band  1. 0-3.0  Hz 

(3) :  Event  not  observed 
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Table  2-10:  Magnitude  Yield  Relations  for  Kazakhstan  from 
6  Events  Recorded  at  the  GERESS  Site  1988-1989 


Yr 

Mo 

Day 

Rem 

“b 

GRF 

Yield 

GRF 

mb 

GERESS 

Yield 

Nuttli 

Yield 

Ringdal 

Yield 

Ekstrom 

1988 

May 

4 

6.0 

150 

6.0 

147 

132 

127 

1988 

Jun 

14 

4.6 

10 

4.4 

4 

7 

1988 

Sep 

14 

JVE 

6.0 

149 

5.9 

108 

108 

119 

1988 

Dec 

17 

5.7 

71 

5.8 

80 

68 

131 

1989 

Jan 

22 

6.0 

120 

6.0 

147 

118 

111 

1989 

Feb 

12 

5.8 

90 

6.0 

147 

63 

88 

The  yields  of  column  Nuttli  were  calculated  using  his  1986  paper  and  GERESS  magnitudes. 
Yields  in  columns  Ringdal  are  from  Ringdal  et  al.  (1992),  and  Ekstiom  from  Ekstrom  and 
Richards  (1994). 


34 


35 


Figure  2-1:  Location  of  the  GERESS  array, 
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Figure  2-3:  Location  of  the  nuclear  test  sites  monitored  at  GERESS  (1988  - 
June  1994). 
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Figure  2-4:  Optimum  GERESS  beams  (filtered  0.5  -  2.5  Hz)  of  events  at 
Nevada  test  site  (1990-1993),  sorted  according  to  magnitude  determined  at 
GERESS  (Table  2-2).  For  the  2  events  in  1988,  single  traces  are  shown.  All 
traces  start  10  seconds  before  the  theoretical  onset  time  calculated  from  the 
published  DOE  parameters  and  the  IASP91  travel  time  tables. 
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Figure  2-6:  Same  as  Fig.  2-4,  except  filter:  1.0  -  3.0  Hz. 


>iTelaiiori  matrix  of  the  eveats  from  Nevada 
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test  site  (for  event  numbers,  see  Table  2-3).  Note  that  the 
coefficient  is  only  0.34  for  event  2. 
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Figure  2-9:  Map  of  event  locations  at  Nevada  test  site  (using  DOE  informa¬ 
tions,  NEIC), 
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Tuamotu  Archipelago  in  1990-1991  (Tables  2-4  and  2-5).  For  the  2  events  in 
1988,  single  traces  are  shown.  All  events  were  sorted  according  to  yield  esti¬ 
mate.  The  trace  of  the  event  on  Nov.  14,  1990  starts  10  seconds  before  the 
measured  onset  time,  the  other  traces  were  aligned  to  its  first  minimum  (the 
start  time  of  each  trace  is  indicated  to  the  right). 


YIELD  ESTIMATE  (DSIR)  [kt] 


Figure  2-11:  Relation  of  mbCPKPcERESs)  versus  estimated  yield  for  events  at 
Tuamotu  Archipelago. 
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Figure  2-12:  Cross-correlation  sequence  of  the  optimum  beam  for  the  event  on 
Nov.  21,  1990  in  the  Tuamotu  Archipelago  with  the  optimum  beams  of  all 
other  events  (1990-1991). 
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Figure  2-13:  Coirelaiion  matrix  of  the  events  from  Tuamotu  Archipelago  (for 
event  numbers,  see  Table  2-5). 
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EVENT  NUMBER  (91.7.5  =  1. 90.11.14  =  14) 
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Figure  2-15:  Optimum  GERESS  beams  (filtered  0.5  -  2.5  Hz)  of  events  at  Lop 
Nor,  sorted  according  to  magnitude  determined  at  GERESS  (Table  2-7).  The 
trace  of  the  event  in  1993  starts  10  seconds  before  the  measured  onset  time,  the 
other  traces  were  aligned  to  its  first  maximum  (the  start  time  of  each  trace  is 
indicated  to  the  right). 


05.00. 52. 57S 


Figure  2-16:  Untiltered  Deams  ot  f  ana  fcf  rrom  me  i^p  iNor  eveni  on  ma; 
21,  1992  (top)  and  the  corresponding  true  ground  motion  simulations  (bottom). 
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Figure  2-17:  Optimum  GERESS  beam  (filtered  0.5  -  2.5  Hz)  of  the  1990  event 
on  Novaya  Zemlya  (Table  2-8).  The  trace  starts  at  15:04:04.0. 
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3.  P-WAVE  ARRAY  POLARIZATION  ANALYSIS  AND 
EFFECTIVE  ANISOTROPY  OF  THE  BRITTLE  CRUST* 


G.  H.R.  Bokelmann 


3.0  ABSTRACT 


P-wave  polarization  constrains  local  anisotropy  in  the  vicinity  of  the  re¬ 
ceivers.  Using  three-component  and  array  data  from  the  regional  GERESS 
array  in  Southeastern  Germany,  we  measure  polarization  a  and  propaga¬ 
tion  vectors  s  for  P-phases  of  120  events.  Angular  deviations  a-s  between 
these  normalized  vectors  often  approach  10®,  rendering  them  easily  measur¬ 
able.  The  effect  of  anisotropy  can  be  distinguished  from  remote  effects,  since 
all  remote  effects  such  as  source  mislocation,  distant  lateral  heterogeneity 
or  distant  anisotropy  affect  polarization  and  propagation  vectors  simiiltane- 
ously.  Averaging  removes  sensitivity  to  near-receiver  heterogeneity,  and  local 
anisotropy  is  left  as  sole  cause  of  the  effects  in  a-s.  This  method  hence  gives 
local  effective  anisotropy  in  the  neax-receiver  crust  averaging  over  a  depth 
interveJ  of  a  wavelength  (w  6km).  We  resolve  strike  and  dip  of  the  symetry 
plane  and  also  two  dimensionless  numbers  t)  and  r  which  give  constraints  on 
four  of  the  elastic  parameters.  The  optimtim  model  (variance  reduction  44%) 
has  symetry  plane  orientation  of  strike  113®  and  dip  49®  to  the  North,  which 
corresponds  closely  to  the  consistently  observed  gneiss  foliation  direction  in 
the  area  (120®,  50®-60°).  Comparing  t]  and  r  with  predictions  from  different 
physical  models  we  find  that  the  data  axe  fit  by  a  gneiss  model  assuming  that 
the  anisotropy  is  dominated  by  the  mica,  if  3%  to  8%  of  the  mica  axe  well- 
aligned.  This  suggests  that  anisotropy  in  the  region  studied  is  dominated  by 
the  effect  of  local  foliation  rather  than  the  regional  stress  field. 

Key  words:  array  seismology,  wave  propagation,  polarization,  anisotropy, 
inversion 

3.1  INTRODUCTION 


In  the  pioneering  study  of  anisotropy  in  the  Earth  Hess  (1964)  showed 
that  azimuthal  travel  time  delays  of  oceanic  P„  can  be  explained  by  aligned 
olivine  crystals  in  the  oceanic  upper  mantle.  For  such  a  phase,  which  travels 
subhorizontally  on  most  of  its  path.  Backus  (1965)  showed  that  azimuthal 
travel  time  variation  is  characterized  by  a  rather  simple  trigonometric  ex¬ 
pression,  given  that  anisotropy  is  weak.  Fitting  of  such  trigonometric  terms 
has  subsequently  been  used  for  teleseismic  P-waves,  f.e.  in  Dziewonski  and 

*  submitted  to  Geophysical  Journal  International 
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Anderson  (1983).  Interpretation  of  such  data,  however,  is  complicated  by  the 
fact  that  lateral  heterogeneity  can  also  cause  azimuthal  variation.  This  trade¬ 
off  is  often  difficult  to  resolve  from  travel  time  data  alone.  The  second  type 
of  ’’classical”  anisotropy  studies  is  shear- wave  splitting  (Ando,  Ishikawa  and 
Yamazaki,  1983;  Silver  and  Chan,  1991),  a  technique  which  uses  the  strong 
effect  of  anisotropy  on  waveforms  of  shear-waves  due  to  different  speed  of 
the  two  orthogonally  polarized  shear-waves.  The  constraint  on  anisotropy  is 
not  unlike  that  of  travel  times.  Both  give  integral  constraints  only,  namely 
on  bulk  anisotropy  along  the  wave  path.  Hence,  for  both  methods  a  tradeoff 
exists  between  size  of  anisotropy  and  extent  of  anisotropic  path.  More  im¬ 
portantly,  an  anisotropic  region  is  difficult  to  localize  using  these  data.  For 
example,  mantle  anisotropy  can  be  inferred  from  shear-wave  splitting  data 
only  if  observed  effects  are  too  large  to  be  explained  by  crustal  amsotropy 
alone.  This  is  somewhat  imsatisfying,  and  methods  for  independent  deter¬ 
mination  of  anisotropy  in  the  crust  under  the  receivers  should  be  sought. 
This  paper,  in  fact,  demonstrates  that  beside  active  experiments  an  entirely 
independent  constraint  on  near-receiver  amsotropy  is  possible;  in  this  case  it 
follows  from  analysis  of  three-component  array  data. 

Anisotropy  studies  using  travel  times  and  shear- wave  splitting  can  be  per¬ 
formed  with  single-station  data,  in  the  later  case  requiring  a  three- component 
instrument.  For  the  present  study  of  P-wave  phases,  we  use  array  and  three- 
component  data  from  the  GERESS  array.  These  array  data  allow  determi¬ 
nation  of  the  propagation  vector  s  =  k/a;  of  plane  waves  crossing  the  array. 
While  this  particular  feature  of  beamforming  is  well-appreciated  for  signal 
enhancement,  event  detection  and  localization  (Beauchamp,  1975),  this  par¬ 
ticular  strength  has  apparently  not  been  used  in  amsotropy  studies  so  far . 

The  propagation  vector,  given  by  s  =  Vt,  may  be  specified  by  azimuth  ds 
and  horizontal  slowness  p  —  or  alternatively  incidence  angle  <j>g  = 

sin~^{pv)  with  the  local  velocity  v.  On  the  other  hand,  three-component 
stations  allow  determination  of  the  polarization  vector  a,  which  gives  the 
dominant  direction  of  linear  particle  motion,  a  can  be  specified  by  azimuth 
Oa  and  incidence  angle  <f>a-  In  the  following  we  restrict  om:  intention  to 
normalized  s  and  a,  namely  to  azimuth  and  incidence  angles.  This  effectively 
ignores  that  the  array  also  gives  estimates  of  the  true  or  apparent  velocity. 
Neither  of  the  two  azimuth  estimates  and  9a  necessarily  correspond  to  the 
azimuth  of  the  receiver-somce  great  circle,  since  heterogeneity,  anisotropy 
and  other  effects  may  cause  deviations.  The  same  is  true  for  the  incidence 
angles.  For  known  sotirce  location  we  get  predictions  {Bp,  <j>p)  for  a  reference 
Earth  model.  Comparing  {Bs,<i)s)  and  {BaAa)  with  the  predictions  gives  rise 
to  so-called  mislocation  diagrams.  For  the  pmpose  of  event  localization  the 
two  data  types  and  correspondingly  their  nuslocation  patterns  are  typically 
treated  the  same,  where  the  azimuth  is  used  for  localization  pmrposes,  while 
the  incidence  angle  serves  for  identification  of  the  phase  for  regional  events 
and  epicentral  distance  in  the  case  of  teleseismic  events. 

Interestingly,  source  location  procedures  have  so  far  disregarded  that  these 
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two  observed  quantities,  namely  propagation  and  polarization  vectors,  axe 
principally  different.  While  they  coincide  in  an  isotropic  Earth,  elastic 
anisotropy  in  the  Earth  would  affect  propagation  and  polarization  vectors 
differently.  For  event  localization,  this  means  that  there  is  a  systematic  bias 
between  propagation  and  polarization  vectors,  which  should  be  removed  by 
taking  local  anisotropy  into  account.  This  systematic  bias,  however,  is  the 
prime  quantity  of  interest  in  this  study,  since  it  offers  us  a  new  constraint 
on  anisotropy  in  the  vicinity  of  the  receivers  (in  fact  exclusively  there).  In 
this  study  we  use  polarization  and  propagation  vectors  exclusively.  Travel 
time  and  amplitude  information  is  not  used,  since  that  might  obscure  the 
near-receiver  crustal  anisotropy  with  information  gathered  somewhere  along 
the  ray  paths. 


3.2  OBSERVATIONS 


3.2.1  GERESS  Array 

The  GERman  Experimental  Seismic  System  GERESS  (Harjes,  1990), 
shown  in  figure  3.1,  is  a  regional  seismological  array  with  apertiire  of  about 
4  km  consisting  of  25  stations  with  1  Hz  vertical  short-period  instruments 
(Geotech  GS-13)  sampled  at  40  Hz,  with  4  of  them  also  equipped  with  hor¬ 
izontal  components  (A2,D1,D4,D7;  triangles).  Additional  three-component 
instruments  at  the  C2  location  are  sampled  at  120  Hz  (Geotech  GS-13)  and 
10  Hz  (Geotech  BB-13).  For  further  descriptions  of  installation  criteria,  per¬ 
formance  and  operation  see  Harjes  (1990)  and  Harjes  et  al.  (1993).  The 
array  is  located  in  Eastern  Bavaria,  Germany,  in  a  hilly  region  on  crystalline 
outcrop  of  the  Bohemian  massive.  The  topographic  variation  across  the  ar¬ 
ray  is  about  200  meters.  A  major  advantage  of  data  from  the  GERESS  array 
is  that  there  are  no  sediments  imder  the  stations,  greatly  improving  the  qual¬ 
ity  particularly  of  the  polarization  data.  Otherwise  P-S  conversions  at  the 
crust-sediment  boimdary  might  disturb  polarizations  substantially. 

3.2.2  Events 

This  study  uses  a  set  of  95  events  recorded  at  GERESS  array  stations  in 
the  time  interval  between  May  1990  and  January  1993  (table  3.1).  These 
events  comprise  the  whole  range  from  regional  to  teleseismic  distance.  The 
emphasis  in  this  study  is  on  first-arriving  phases.  Only  the  later-arriving 
phase  PcP  is  included  in  three  cases.  Table  3.1  is  part  of  a  larger  set  of 
events.  For  the  events  not  included  in  the  table  either  the  polarization  or 
propagation  vector  had  to  be  discarded,  which  was  the  case  particularly 
often  early  in  the  experiment.  In  general,  events  were  discarded  if  data 
from  less  than  3  three-component  stations  or  9  array  stations  were  available. 
Out  of  112  events  with  acceptable  polarization  data,  81  had  data  for  all  6 
three-component  instruments  present.  Out  of  108  events  with  acceptable 
propagation  data,  88  had  data  from  at  least  20  array  stations. 

Next  we  discuss  procedures  of  extracting  backazimuth  and  incidence  angle 
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information  from  three-component  arid  array  data. 

3.2.3  Determination  of  Polarization  Directions 
For  three-component  data  u-^  =  (uiY  with  components  j  =  1,2,3,  the 
covariance  matrix  at  zero-lag  is  computed  as 

Sjk  =  (1) 


where  N  is  the  number  of  samples  in  the  time  window.  Obviously,  S  is 
symmetric,  so  that  its  eigenvalues  A-^  are  real  and  nonnegative.  We  choose 
a  set  of  3  orthonormal  eigenvectors  sj.  \^aP  (no  summation)  then  gives  the 
axis  of  the  polarization  ellipsoid  (see  f.e.  Kanasewich,  1981).  K  the  particle 
motion  is  linearly  polarized,  then  A^  =  A®  =  0.  For  elliptical  particle  motion 
(within  a  polarization  plane),  we  have  A^  =  0.  In  general,  however,  the 
particle  motion  is  ellipsoidal  due  to  presence  of  different  types  of  noise.  If  we 
take  components  (1,2,3)  as  (Vertical  up.  North,  East),  the  incidence  angle 
is  =  cos~^|ui^l  and  the  backazimuth  is 


^itan  for  <  0 

(  TT -f  tara“^(^|r)  for  ai^  >  0 

A  simple  check  of  the  quality  of  the  measurement  is  given  by  the  rectilinearity 


2A1 


(2) 


which  is  equal  to  1  for  linear  polarization  and  0  in  the  worst  case  of  ’’spheri¬ 
cal”  polarization,  when  all  eigenvalues  have  the  same  size. 

Figure  3.2  shows  a  seismogram  example  demonstrating  a  few  characteris¬ 
tics  of  the  eigenvector  analysis.  For  this  data  example  (year  1992  day  48  horn 
19)  with  an  epicentral  distance  of  385km  and  backazimuth  of  235°  a  number 
of  regional  phases  are  visible.  The  main  phases  for  this  distance  range  are 
Pn,  Pg:  5„  and  Sg.  They  are  visible  in  the  three-component  seismogram 
on  the  top.  The  lower  traces  give  polarization  attributes  computed  using 
a  ghding  window  with  a  length  of  2  seconds:  The  root-mean-square  ampli¬ 
tude  of  the  three-component  seismogram  or  ’’three-component  size”  (fourth 

trace)  A^'  =  Sjj  shows  Pg,  Sn  and  Sg  clearly.  The  first- arriving 

shows  up  with  somewhat  smaller  three-component  size,  but  with  the  high¬ 
est  rectilinearity  (2),  with  close- to-ideal  linear  polarization.  A  line  gives  the 
critical  rectilinearity  level  of  0.9.  Only  Pn  and  Pg  exceed  this  level.  Since 
the  eigenvalues  are  related  to  energy,  error  bars  can  be  obtained  from  the 
F-test  (see  Bokelmann,  1992).  Sg  shows  up  with  higher  rectilinearity  than 
other  S  phases  and  also  substantial  three-component  size.  The  horizontal- 

vertical  ratio  (H/V)  and  similarly  the  radial- transverse  ratio  (R/T) 
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computed  from  the  elements  of  the  covariance  matrix  are  also  indicative  of 
the  phase  character.  The  values  of  H/V=l  and  R/T=l  are  shown  by  lines 
indicating  polarization  angles  near  45°.  For  the  S-phase  train  H/V  shows 
very  clear  horizontally  polarized  phases  Sn,  Sg  and  Lg,  which  also  show 
STpall  values  of  R/T.  This  radial-transverse  ratio  for  P  phases  shows  clearly 
radially  polarized  energy.  In  summary,  polarization  is  an  important  tool  for 
phase  identification.  In  the  following  we  want  to  use  polarization  data  for 
the  set  of  P  phases.  Figure  3.2  shows  that  polarization  attributes  for  P- 
phases  are  relatively  slowly  varying.  In  fact,  they  are  also  insensitive  to  the 
particular  choice  of  window  length.  On  the  6  three-component  stations  we 
observe  for  P„  a  backazimuth  of  da  =  240.4°  ±  5.5°  and  an  incidence  angle 
of  =  44.9°  ±  4.2°.  These  angle  estimates  correspond  closely  to  values 
extracted  manually  from  hodograms  and  are  near  predicted  values  from  the 
IASP91  Earth  model  and  source  locations  from  table  3.1  {6p  =  235.17°  and 

=  40.02°). 

Typically  we  have  data  from  six  three-component  stations  available.  We 
computed  standard  errors  simply  from  the  empirical  distribution  of  those 
values.  The  average  uncertainty  for  all  phases  in  the  study  is  3.5°  for  the 
incidence  angle  and  4.9°  for  the  backazimuth,  where  we  used  only  events  with 
high  Signal-to-noise  ratio. 

A  slightly  different  approach  was  given  by  Jurkewicz  (1988),  who  proposed 
averaging  covariance  matrices  of  individual  three-component  stations  before 
performing  the  eigenvector  analysis  to  stabilize  the  estimates. 

Polarization  analysis  lends  itself  to  very  efficient  semi-automatic  process¬ 
ing.  Caution  is  required,  however,  for  phases  closely  spaced  in  time,  particu¬ 
larly  if  their  amplitudes  differ  substantially.  For  regional  phases  this  case  is 
not  imcommon.  Window  lengths  in  this  study  are  allowed  to  vary  depending 
on  data  type  and  quality.  The  median  window  length  is  near  3.5  seconds. 
Similarly,  comer  frequencies  of  a  zero-phase  bandpass  filter  were  allowed  to 
vary  within  0.6  to  4  Hz  to  adapt  to  the  different  wave  types. 

3.2.4  The  Mislocation  Hemisphere 

Figure  3.3  shows  part  of  the  lower  hemisphere  (0°  to  50°)  of  polarization 
anomalies,  namely  deviations  of  observed  from  predicted  polarization,  where 
predicted  again  refers  to  values  computed  from  IASP91  model  with  source 
locations  from  table  3.1.  Here  the  free-surface  effect  is  taken  into  accotmt, 
as  wiU  be  discussed  below.  In  such  a  display,  the  center  corresponds  to 
vertical  incidence.  The  perimeter  gives  50°  incidence  with  North  at  the 
top.  We  choose  this  type  of  diagram  to  fa.cilitate  comparison  with  predicted 
polarization  patterns  from  anisotropic  models. 

Multiple  nuclear  events  from  two  source  regions  (Nevada,  Tuamotu)  give 
rather  consistent  deviations  of  polarization  to  Northern  directions.  Also 
other  geographic  regions  with  several  nearby  events  give  rise  to  consistent 
polarization  effects.  For  steeply  incident  phases  from  the  North  there  is 
a  general  deviation  to  more  Northern  directions,  while  the  Southern  part 
of  figure  3.3  shows  tangential  deviations,  apparently  varying  rapidly  over 
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the  lower  hemisphere.  For  example  polarization  anomalies  from  Northern 
Italy  and  Switzerland  show  opposite  tangential  effects.  However,  based  on 
polarization  data  alone,  this  behaviour  can  not  be  interpreted  in  terms  of 
heterogeneity  or  anisotropy.  Consistently  occuring  source  mislocation  effects 
for  these  regional  distances  might  also  cause  such  variation.  While  this  may 
be  a  problem  of  either  polarization  vector  data  or  propagation  vector  data 
tahen  by  itself,  the  difference  is  insensitive  to  source  location  errors  and  nei¬ 
ther  to  distant  lateral  heterogeneity.  Instrument  calibrations  were  performed 
daily.  Substantial  deviations  from  the  nominal  response  were  rare  and  the 
correction  had  a  rather  small  effect  on  the  lower  hemisphere  representation. 

3.2.5  Free  Surface 

Conversion  at  the  free  surface  perturbs  the  incidence  angle  of  the  polar¬ 
ization  direction.  The  apparent  incidence  angle  a  for  a  P-wave  is  given 

by  _ 


$'a  =  tan 


2sin{^a)\/A'^  — 

—  2sin2(#a) 


(3) 


with  the  velocity  ratio  A  —  Vp/vg  (Aid  and  Richards,  1981).  The  dependence 
of  on  $  is  shown  in  figure  3.4  for  a  number  of  Poisson’s  ratios  «,  which 
are  related  to  the  velocity  ratio  A  by  «;  =  .  The  three  values  are 

K  =  0.2  (A  =  1.63;  dotted  line)  ,k  =  0.25  (A  =  1.73;  sohd  fine)  and  k  =  0.3 
(A  =  1.87;  broken  fine).  Figvure  3.4  shows  that  the  free-surface  effect  causes 
incidence  angle  deviations  of  several  degrees.  The  free-surface  correction 
is  particularly  important  for  large  incidence  angles.  Variation  of  k  from 
K  =  0.25  to  0.2  or  0.3  produces  offsets  within  5®  in  the  incidence  range 
of  our  study  (4°  to  40°).  The  deviation  between  the  curves  is  within  a 
few  degrees.  For  the  free-siirface  correction  (figure  3.3)  we  used  a  crustal 
Poisson  ratio  of  /c  =  0.25.  This  corresponds  closely  to  the  upper  crustal 
mean  value  observed  in  the  KTB  drill  hole,  which  is  located  about  150km 
northwest  (figure  3.1a)  in  a  similar  geological  setting  (figure  3.1).  Without 
this  correction,  the  lower  hemisphere  has  a  distinctly  different  appearance 
(figure  3.5)  with  a  clear  radial  component.  The  mean  of  this  radial  component 
is  5.1°  oriented  outwards.  Clearly,  the  free-surface  correction  is  an  important 
and  observable  effect,  which  must  be  taken  into  account.  On  the  other  hand, 
we  have  seen  that  there  are  clear  effects  in  figure  3.3  beyond  the  simple  effect 
of  crustal  Poisson  ratio  alone.  For  the  Southern  part  of  figure  3.3  the  radial 
effect  almost  vanishes. 

3.2.6  Array  Meeisurements 

Propagation  directions  (^^,$3)  can  be  extracted  from  array  data  either 
from  travel  times  or  from  fk-analysis.  Characteristics  of  this  problem  have 
been  studied  by  Bokelmann  (1993),  where  it  was  shown  that  small-array  data 
like  the  ones  in  this  paper  can  produce  propagation  vector  data  with  reason¬ 
able  imcertainty  given  that  topography  (figure  3.1)  is  taken  into  accotmt.  We 
obtain  the  mislocation  pattern  of  figure  3.6.  Conversion  from  slowness  p  to 
incidence  angle  =  sin~^(pv)  was  done  using  the  local  velocity  estimate  of 
v=5.2  km/s  obtained  from  3D  array  analysis  (Bokelmann,  1993).  In  that  pa- 


58 


per,  the  reader  also  finds  a  discussion  on  the  consistency  of  propagation  vector 
estimates  from  examples  of  Nevada  and  Tuamotu  nuclear  events.  These  data 
are  included  in  this  data  set.  In  figure  3.6  they  show  comparatively  small 
deviation  dociimenting  consistency  for  the  array  data.  Nevertheless  these 
deviations  showed  up  as  statistically  significant  (Bokelmann,  1993).  There 
are  even  larger  anomalies  in  these  data  though.  The  Northern  deviation  for 
steeply  incident  phases  from  the  North  which  was  observed  for  polarization 
data  (figure  3.3)  is  not  apparent  in  the  array  data.  Interestingly  the  pattern 
for  Yugoslavia  and  Romania  are  very  similar  to  observed  polarization  anoma¬ 
lies  in  figure  3.3  as  well  as  the  rapid  variation  of  Switzerland  and  Northern 
Italy.  This  suggests  that  either  lateral  heterogeneity  or  source  mislocation  is 
affecting  propagation  and  polarization  directions. 

For  an  epicentral  distance  of  about  400km  the  mean  tangential  deviation 
between  Switzerland  and  Northern  Italian  events  would  require  a  relative 
source  mislocation  of  about  70km  occuring  consistently  in  the  PDE  solutions 
(Preliminmy  Determination  of  Epicenters)  between  the  two  source  regions. 
While  this  is  not  impossible,  lateral  heterogeneity,  the  Alpine  root,  might  also 
cause  the  azimuthal  deviation,  or  respectively,  the  relative  source  mislocation. 
Both  would  simultaneously  perturb  propagation  and  polarization  directions. 

Nevertheless,  the  procedure  of  relative  analysis  of  polarization  and  prop¬ 
agation  data  gives  data  which  are  independent  of  the  source  location  and 
distant  lateral  heterogeneity.  The  relative  data  (polarization  -  propagation 
vectors)  from  the  94  events  are  shown  in  figure  3.7.  However,  Nevada,  Tu¬ 
amotu,  and  other  regions  still  show  up  with  consistent  effects.  Naturally,  the 
random  scatter  is  larger  for  these  differences  than  for  the  individual  data.  For 
two  events  the  difference  data  were  clearly  dominated  by  noise  and  had  to 
be  removed.  It  may  seem  surprising  that  source  mislocation  and  lateral  dis¬ 
tant  heterogeneity  don’t  affect  the  relative  data.  In  fact,  the  method  doesn’t 
require  knowledge  of  source  locations  at  all.  The  propagation  vector  is  used 
to  determine  the  location  in  the  lower  hemisphere.  Measmred  from  this  di¬ 
rection,  the  polarization  vector  deviation  indicated  presence  of  anisotropy. 
This  allows  studying  anisotropy  imder  the  receivers.  Lateron,  we  use  the 
data  of  figme  3.7  in  an  inversion  for  anisotropy.  But  first,  we  review  wave 
propagation  theory  in  anisotropic  media. 

3.3  REVIEW  OF  PROPAGATION  IN  ANISOTROPIC  MEDIA 


Here  we  wish  to  merely  state  a  few  features  of  particular  importance  for 
this  study.  Comprehensive  treatment  of  the  topic  has  been  given  f.e.  by 
Crampin  (1981). 

In  a  general  anisotropic  medium,  the  solution  for  a  propagating  plane  wave 
front  follows  from  the  Ansatz 

u(r,  t)  =  a(r)/(t  -  nr/ c)  (4) 

where  a(r)  describes  amplitude  and  polarization  of  the  wave,  while  the  time- 
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dependence  is  given  by  /()  depending  on  the  phase  velocity  c.  From  the 
equation  of  motion,  we  get 

muai  =  c^Oi  (5) 

with  the  Christoffel  matrix 

mu  =  CijkinjUklp  (6) 

(6)  has  the  form  of  an  eigenvalue  equation  with  the  square  of  the  phase  ve¬ 
locity  appearing  as  eigenvalue.  Hence  polarization  vectors  a  serve  as  eigen¬ 
vectors.  Due  to  the  symmetry  of  the  elasticity  tensor  Ciju,  we  can  always 
choose  a  set  of  3  mutually  perpendicular  eigenvectors.  Below  we  will  give 
examples  of  Ciju- 

The  tj'pical  way  to  solve  (6)  is  to  specify  a  propagation  direction  n  and 
solve  for  eigenvalue  and  eigenvector  for  each  of  the  three  polarizations  (in 
the  nondegenerate  case).  It  is  clear  that  in  general  propagation  and  polar¬ 
ization  vectors  don’t  coincide  for  waves  propagating  in  anisotropic  media. 
For  isotropic  media  that  was  the  case.  Therefore  we  can  use  such  deviations 
to  determine  effective  medium  anisotropy. 

For  the  case  of  P-waves,  we  will  discuss  this  discrepancy  more  closely 
below,  since  both  of  these  directions  can  be  observed  using  arrays  and  three- 
component  stations. 

The  wavefront  propagates  in  the  normal  direction  with  the  phase  velocity 
c  =  |Vf|“^.  Vf  is  the  propagation  vector,  which  can  be  estimated  from  data 
of  small  arrays  (see  also  Helbig,  1958). 

The  polarization  vector  a  gives  the  direction  of  local  particle  motion.  It 
is  measured  directly  by  polarization  analysis.  Neither  of  these  two  vector 
should  be  confused  with  the  group  velocity  vector,  which  is  the  direction, 
in  which  energy  propagates  (Musgrave,  1970).  Crampin,  Stephen  and  Mc- 
Gonigle  (1982)  showed  that  for  weak  anisotropy  the  difference  between  po¬ 
larization  and  group  velocity  directions  are  small.  This  is  intuitively  clear, 
since  the  energy  should  propagate  in  a  direction  near  the  particle  motion. 

While  the  angular  deviation  between  phase  and  group  velocity  vectors 
received  quite  some  interest  in  the  literatme  (Crampin,  1981;  Crampin  et 
al.,  1982),  relation  between  phase  velocity  and  polarization  vectors  has  been 
discussed  for  model  crystals  (Crampin,  1978,  1981;  Crampin  et  al.,  1982), 
but  due  to  a  lack  of  appropriate  data,  there  are  only  few  applications  to  real 
data  so  far  (White  et  al.,  1983;  Li,  Leary  and  Aki,  1987;  de  Parscau,  1991a). 
However,  deviations  of  polarization  directions  from  propagation  directions 
can  be  substantial  as  we  will  see  below.  In  the  following  we  want  to  illustrate 
this  effect  and  lateron  exploit  it  in  an  inversion  for  effectively  anisotropic 
structure.  It  is  important  to  note  that  the  polarization  depends  only  on 
local  anisotropy.  This  is  illustrated  in  figure  3.8. 

In  the  context  of  anisotropy,  P-wave  polarizations  have  been  studied  for 
VSP  borehole  data  by  de  Parscau  (1991a)  and  Li,  Leary  and  Aki  (1987). 
However,  with  a  ID  array  in  the  borehole,  there  is  no  control  over  off-azimuth 
arrivals. 
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We  have  seen  that  our  method  of  relative  analysis  is  insensitive  against 
distant  heterogeneity.  But  how  about  near-receiver  heterogeneity?  Inspec¬ 
tion  of  polarization  hemispheres  for  individual  three-component  instruments 
showed  some  station-dependent  scatter  which  is  partly  due  to  deviations 
of  the  instrument  response  from  nominzd  values.  Small-scale  heterogeneity 
near  the  receivers  may  also  cause  waveform  distortions.  Hence,  polariza¬ 
tion  data  from  single  three-component  stations  have  to  be  regarded  with 
some  caution.  Such  data  carry  information  about  the  local  conditions  imder 
the  station.  For  our  piirpose  we  average  over  all  three-component  stations 
of  the  array,  thereby  obtaining  average  polarization  attributes,  where  array 
stations  are  well- distributed  over  an  aperture  of  4  km.  The  effect  of  near¬ 
receiver  small-scale  heterogeneity  is  expected  to  be  nearly  uncorrelated  from 
station  to  station  and  we  expect  the  averaging  procedure  to  be  effective  in 
averaging  out  most  of  its  effect  on  the  difference  pattern.  The  effect  from 
larger-scale  heterogeneity  on  the  other  hand,  is  likely  correlated  from  station 
to  station.  We  obtain  a  net  effect  in  the  polarization  vectors,  but  simi¬ 
larly  in  the  propagation  vectors.  The  difference  should  therefore  again  be 
approximately  unaffected.  This  argument  suggests  that  although  scatter  is 
im'ariably  present  in  difference  patterns  hke  figure  3.7,  the  scatter  is  domi¬ 
nated  by  random  fluctuation  and  not  by  biases  from  discrete  local  scatterers. 
A  very  important  feature  of  the  GERES S  array  is  the  absence  of  a  seismo- 
logically  significant  near-surface  low- velocity  layer,  what  can  be  shown  from 
inspection  of  the  train  of  P-S  conversions. 


3.4  UPPER-CRUSTAL  ANISOTROPY 


Anisotropy  in  the  vicinity  of  the  receivers,  however,  does  affect  the  relative 
data  in  figure  3.7.  Which  depth  range  is  it  sensitive  to?  We  have  seen  that 
polarization  is  a  local  property.  For  seismological  phases  of  finite  wavelength 
this  means  that  the  depth  range  equivalent  to  about  a  wavelength  should 
affect  the  observed  polarization.  For  the  typical  dominant  frequencies  in  this 
study  we  have  wavelengths  of  about  6  km.  Hence,  we  expect  to  estimate 
a  spatial  average  of  anisotropic  properties  of  the  top  few  kilometers  in  the 
upper  crust.  Since  we  average  over  all  three-component  stations  in  the  array, 
we  assume  that  discrete  heterogeneity  within  the  array  can  be  ignored  and 
allow  for  an  anisotropic  but  homogeneous  model. 

In  this  context,  there  are  several  phenomena  which  may  cause  effective 
anisotropy.  Intrinsic  anisotropy  of  minerals,  as  well  as  small-scale  structure 
possessing  certain  symetry  lead  to  directional  dependence  of  macroscopic 
elastic  parameters  and  hence  propagation  velocity.  Sedimentary  layering, 
with  a  potentially  prominent  effect,  is  not  present  here.  However,  cracks, 
fractures  and  rock  foliation  (Babuska  and  Cara,  1991),  if  certain  orientations 
prevail,  can  give  rise  to  substantial  effects.  Since  our  data  are  sensitive  to 
a  spatial  average  over  several  kilometers,  we  expect  macroscopic  anisotropic 
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behaviour  either  if  the  anisotropy  is  locally  very  strong  or  if  preferred  orien¬ 
tation  is  spatially  coherent  over  a  depth  range  of  several  kilometers  within 
the  topmost  part  of  the  upper  crust.  The  latter  may  for  example  be  the  case 
for  cracks  caused  by  the  regional  stress  field.  Note  that  our  P-wave  polariza¬ 
tion  data  are  sensitive  only  to  the  depth  range  where  brittle  deformation  is 
expected. 

General  anisotropic  model/hexagonal  case 

Since  the  cause  of  anisotropy  is  not  known  beforehand,  it  is  preferable 
to  find  elastic  constants  imder  no  restricting  assumptions.  In  the  follow¬ 
ing  we  will  study  the  type  of  constraints  P-wave  polarization  data  pose  on 
elastic  constants.  We  wiU  restrict  our  attention  to  the  case  of  hexagonal 
symetry,  since  this  is  a  viable  model  for  most  of  the  possible  causes  of  ef¬ 
fective  anisotropy  mentioned  above.  An  exception  may  be  more  complex 
intrinsic  anisotropy  of  individual  crystals,  which  is  not  considered  here.  In 
the  case  of  a  hexagonal  or  transversely  isotropic  symetry,  the  elastic  constant 
matrix  Cijki  has  the  simple  shape 


c-ijki  = 


/  Cii 

Cii  2cgg 

Cl3 

\ 

Cl  a  —  2c66 

Cll 

Cl3 

Cl  3 

Cl3 

C33 

C44 

^44 

\ 

C66  / 

(7) 


with  5  unknowns,  assuming  that  the  symetry  axis  is  in  xg-direction.  Given 
this  particular  symetry,  velocity  v  and  polarization  direction  of  P- waves  de¬ 
pend  only  on  the  angle  x  between  propagation  direction  and  symetry  axis 

3.4.1  Modelling  of  P- Velocity  and  Polarization  Patterns 

White  (1983)  and  de  Parscau  (1991b)  give  exact  expressions  for  velocity 


PV^  =  -(Cii  -f  C44)sm^X  +  (C33  +  C44)C0S^X+ 


(  (cil  -  C44)5m^X  -  (C33  —  C44,)cOS^xY  +  4(ci3  -|-  €44)^ cos^ x^^^^ X 
and  polarization  angle  ^  measured  from  the  symetry  axis 

(ci3  -I-  C44)cosx5mx 


1/2 


tan^  = 


pv^(x)  -  (ciisin^x  +  C44Cos2^) 


(8) 


(9) 


(8)  and  (9)  depend  only  on  the  fom  elastic  constants  cn,  C13,  C33  and  C44. 
The  last  constant  cee  in  (7)  can  be  determined  only  from  quasi-transversely 
polarized  shear  waves,  de  Parscau  (1991b)  showed  that  (9)  can  be  simplified 
to 


tan2^  = 


T]sin2x 


1  —  (1  -k  T)sin^x 


(10) 
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using  the  two  parameters 


Cl3  +  C44 

77  = 

(11) 

C33  —  C44 

Cii  *  C44 
r  = - 

(12) 

C33  —  C44 

This  shows  that  we  will  i.g.  not  be  able  to  independently  determine  the 
four  elastic  constants  Cn,  C13,  C33  and  C44  from  P-polarization  data.  Instead, 
constraints  on  these  constants  must  be  derived  from  the  77  and  r  using  the 
nonlinear  relations  (11)  and  (12).  In  addition  to  77  and  r,  the  orientation 
of  the  symetry  axis  is  imknown.  Therefore  we  have  2  additional  unknowns 
representing  strike  and  dip  of  the  symetry  axis.  Both  r  and  77  can  be  seen 
as  measures  of  anisotropy.  Since  cn  and  C33  give  pv^  parallel  and  normal  to 
the  symetry  planes,  r  gives  a  measure  of  the  velocity  deviation  between  these 
special  directions.  77  is  more  difficult  to  interpret,  but  it  can  be  approximately 
related  to  the  P- velocity  anisotropy  at  small  incidence  angles  nearly  normal 
to  the  symetry  planes  (de  Parscau,  1991b).  For  an  isotropic  medium  we  have 
77  =  T  =  1. 

A  number  of  criteria  on  the  elastic  coefficients  in  (7)  can  be  derived  from 
requiring  stability,  namely  the  condition  that  deformation  is  associated  with 
nonnegative  internal  energy.  The  requirement  is  that  (7)  is  positive  semidef- 
inite  (Backus,  1962).  It  follows  that  all  diagonal  elements  axe  nonnegative. 
Furthermore,  cn  >  ces,  033(011  -  cee)  >  ci3^  In  our  context,  bounds  on  elas¬ 
tic  coefficients  are  important  for  ffiading  boimds  on  the  parameters  77  and  r 
and  for  rejecting  physically  impossible  models.  Postma’s  (1955)  inequality 

(cii  —  C44)(C33  —  C44)  >  (ci3  -f  044)^  (13) 

shows  that  r  >  0.  Since  r  gives  a  measme  of  the  size  of  P-wave  anisotropy,  an 
upper  bound  may  be  given  from  an  estimate  of  maximum  possible  amsotropy. 
The  sign  of  cis  -f  C44,  the  nominator  of  77  is  foimd  to  have  substantial  impor¬ 
tance  for  the  polarization  behaviour.  Cases  with  C13  4-  C44  <  0  are  physically 
possible,  but  they  are  foimd  only  is  very  carefully  planned  laboratory  ex¬ 
periments  (Helbig  and  Schoenberg,  1988).  In  these  cases,  the  polarization 
angle  is  found  to  rotate  in  a  sense  opposite  to  the  propagation  direction,  in 
contrast  to  the  normal  case  with  rotation  in  the  same  sense.  For  real  Earth 
studies  it  is  reasonable  to  assume  C13  -f  C44  >  0,  which  gives  us  a  lower  bound 
77  >  0.  An  upper  boimd  on  77  can  be  found  by  applying  a  requirement  that 
C11/C44  and  C33/C44  be  within  a  certain  range. 

Now  we  wiU  illustrate  polarization  behaviour  for  two  special  cases,  namely 
(77,  t)=  (1.43,  2.13)  and  (0.82, 1.00).  These  cases  correspond  to  elastic  param¬ 
eters  for  models  HCDl  and  HCSl  (Crampin,  1984),  which  give  the  effective 
anisotropy  associated  with  a  model  of  oriented  cracks  (Hudson,  1981)  in  the 
liirdting  cases  of  fluid-saturated  and  dry  conditions.  To  derive  this  model 
of  oriented  cracks,  Hudson  (1980)  used  mean  fleld  theory  (Keller,  1964)  and 


63 


derived  expressions  valid  to  second  order.  In  his  perturbative  treatment  the 
stiffness  matrix  can  be  represented  as 

Cijkl  —  Cijkl  "h  Cijicl  d"  Cijkl 

for  long  wavelengths,  where  the  zero-order  term  gives  the  isotropic  reference 
model  specified  by  the  Lame  parameters  A  and  // 


=  XSijSki  +  fi{SikSji  +  SiiSjk)  (15) 


The  first  and  second-order  perturbation  terms  can  be  interpreted  as  first- 
order  effect  of  the  cracks  and  crack-crack  interaction.  For  the  crack  density 
e  =  Na^  /V  (N/V=number  of  cracks  per  volume,  a=mean  crack  radius)  they 
axe 


Cijkl 
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/  A2  A2 

A2  A2 

£  A(A  -j-  2//)  A(A  -f-  2//) 
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A(A  -|-  2^)  \ 
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(17) 

with  q  =  15(-)^  d-  28(-)  d-  28  and  x  =  2^^^.  Here  the  cracks  are  aligned 
perpendicular  to  X3 .  f)  is  the  diagonal  matrix  with  trace 


trace(D)  =  (uii,uii,wii,i^33,  ^33,6) 


(18) 


Ukm  depends  on  conditions  on  the  crack  face.  The  dependence  of  Cijki  on  the 
crack  parameters  is  only  through  D  and  is  fiirther  given  by 


4  (Ad- 2//) 

(19) 

and 

16  (Ad-2/i)  ^ 

3  + 

(20) 

with  k  and  m  given  as 

Trdfi  A  d-  A^) 

(21) 
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and 


4/i^  (A  +  2//) 
^  Trdfi  (3A  +  4/x) 


(22) 


The  e-reinge  of  interest  in  this  study  is  between  0  and  0.1  .  For  such  low  val¬ 
ues  of  e,  the  second-order  approximation  is  considered  uncritical  (Crampin, 
1981).  d  is  the  (large)  aspect  ratio  of  the  cracks.  A'  and  n'  are  the  Lame- 
parameters  within  the  cracks.  The  two  extreme  cases  with  partictdar  physical 
significance  are  1)  /i'  =  A'  =  0  and  2)  /x'  =  0  with  nonzero  A'. 

The  first  case  models  dry  cracks  (HCDl),  with  free  surfaces  at  the  crack 
faces.  In  this  case  k  and  m  vanish.  The  polarization  pattern  is  deter¬ 
mined  solely  by  (77,  r)  =  (1.43,  2.13).  P-wave  phase  velocities  and  polar¬ 
izations  for  this  model  of  vertical  East- West  oriented  cracks  are  shown  in 
figure  3.9  in  lower  hemisphere  representations,  where  for  display  p-urposes 
the  symetry  axis  is  shown  pointing  upwards  in  the  paper  plane.  We  choose 
the  imconstrained  absolute  levels  of  cn,  C44  and  Ces  to  conform  to  model 
HCDl  (backgroimd  velocity  up  =  b.Skm/sec,  vs  =  ZM9km/sec,  density 
p  —  2.6kmf  sec). 

This  choice  affects  only  the  P-velocity  pattern,  but  not  the  polarization 
pattern,  which  depends  only  on  77  and  r.  Propagation  in  the  plane  of  aligned 
cracks  gives  high  velocity  and  the  P-wave  anisotropy  for  this  ideal  model  is 
24%.  The  right-hand  side  shows  polarization  anomalies,  namely  deviations 
of  the  polarization  from  propagation  directions.  In  general  the  polarization 
direction  deviates  into  the  direction  to  larger  phase  velocity,  diverging  at  the 
symmetry  axis  and  converging  to  the  symmetry  plane. 

For  the  second  case  (model  HCSl),  which  models  saturated  cracks,  we 
obtain  figure  3.10.  The  parameters  are  the  same  as  in  figure  3.9,  but  with 
aspect  ratio  <i=0.0001  and  A' =2.25  GPa.  The  polarization  pattern  is  then 
given  by  77=0.82  and  r=l.  Note  that  this  velocity  pattern  differs  substan¬ 
tially  from  the  dry-crack  model:  The  symmetry  axis  has  also  fast  velocity. 
This  gives  rise  to  r=l.  Clearly,  r  is  very  sensitive  in  distinguishing  dry 
crack  conditions  from  saturated  ones.  The  difference  between  minimum  and 
maximiun  velocity  with  3.5%  anisotropy  is  substantially  smaller  than  for  the 
dry-crack  model.  The  maximum  velocity  achieved  is  somewhat  higher  than 
for  the  dry-crack  model  though.  Polarization  anomalies  are  smaller  but  give 
complexer  pattern. 


3.5  MODEL  FITTING 


Observed  polarization  anomalies  can  be  related  to  the  parameters  77  and 
T  describing  the  anisotropy.  The  inverse  problem  requires  solving  the  eigen¬ 
value  problem  (6).  In  this  paper,  this  is  done  by  nonlinear  parameter  fitting. 
As  shown  above,  parameters  in  the  inversion  are  77  and  r,  but  also  two  angles 
6  and  $  describing  azimuth  and  dip  of  the  symmetry  axis.  In  general,  we 
will  not  be  able  to  determine  the  4  elastic  constants  independently.  Instead 
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we  may  obtain  two  nonlinear  constraints  on  these  parameters  from  rj  and 
r,  which  we  may  supplement  with  a  priori  information.  In  our  case,  such 
information  is  given  by  reasonable  assumptions  about  the  reference  velocity. 
In  the  following  we  will  use  reference  velocities  used  in  the  above  examples 
HCDl  and  HCSl.  This  enables  us  to  also  display  velocity  patterns  for  the 
inversion  results. 


For  each  event  we  have  an  observation  of  polarization  direction 


and 


propagation  direction  compare  these  angular  differences 

^  with  predictions  from  models  depending  on  the  pa- 


for  the  i-th  event. 


rameters  (  6,  r],  r  ),  which  we  call  ( 

We  obtain  the  best  fitting  model  by  inspecting  the  misfit  on  the  sphere 


1  ^  r  / 

S(6,  77,  r)  =  ^  f  cos(A$,-  —  A6i)cos(A6i  —  A$/)  j 


(23) 


which  is  the  analogon  of  the  standard  sum-of-squares  misfit  criterion  on  the 
sphere. 

The  models  HCDl  and  HCSl  have  shown  that  the  polarization  pattern 
may  vary  considerably,  depending  on  the  choice  of  77  and  r.  We  start  the 
discussion  of  inversion  results  by  studying  the  dependence  of  the  misfit  on 
symetry  axis  orientation  for  these  two  special  cases.  77  and  r  are  fixed  to  (1.43, 
2.13)  in  the  first  case  of  dry  cracks  and  (0.82,  1.)  in  the  case  of  saturated 
cracks.  We  are  left  with  fitting  6  and  #.  The  nonlinear  search  may  use  a 
rotation  of  the  elastic  constant  matrix 


Cijkl  —  Vip^jg^kr^ls^pgrs  (^^) 

with  the  rotation  matrix  77.  While  this  allows  for  all  types  of  symetry,  we 
may  more  efficiently  use  (9)  in  oxir  Ccise  of  hexagonal  symetry.  Results  for  the 
dry-crack  case  are  given  in  figmre  3.11  showing  the  achieved  misfit  S  given  in 
degree^  by  (23)  depending  on  the  orientation  of  the  symmetry  axis,  namely 
its  location  in  the  lower  hemisphere.  Axes  orientations  with  small  misfit  are 
near  the  circumference.  The  smallest  achieved  misfit  is  135  degree^,  which  is 
substantially  larger  than  the  starting  misfit  of  68.8751  degree^.  The  simple 
dry-crack  model  is  clearly  incompatible  with  the  data.  On  the  other  hand, 
the  saturated-crack  model  (figure  3.12)  gives  reasonable  variance  reduction  of 
about  41%.  Figure  3.12a  shows  that  small  misfits  are  achieved  for  S3rmmetry 
axis  locations  in  two  regions  in  the  Southwest  and  Northeast.  Only  the 
Southwestern  region  gives  rise  to  statistically  significant  solutions.  The  lowest 
contour  shows  the  confidence  region  where  the  critical  misfit  level  is  computed 
from  the  f-distribution  as 

Scrit  =  -SoptU  +  jy  _  ^fk,N-k{^  “  «)]  (25) 
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(Jenkins  and  Watts,  1969),  where  Sopt  is  the  minimxmi  misfit  and  k=4  the 
number  of  parameters.  Higher  contours  give  2  and  4  times  the  separation  of 
the  critical  level  (1  —  a  =  0.95%  confidence  level)  from  the  minimum  level. 
The  optimmn  model  has  a  symmetry  axis  orientation  of  azimuth  6  =  208° 
and  dip  $  =  36°.  The  symmetry  plane  is  oriented  normal  to  the  symmetry 
axis,  and  has  strike  118°  and  dip  54°  to  the  North.  The  optimum  model 
corresponding  to  the  global  minimum  in  the  Southwest  is  shown  in  figure 
3.12b  by  its  velocity  and  polarization  pattern.  To  be  directly  comparable  with 
the  data,  only  the  partial  lower  hemisphere  from  0°to  55°is  shown.  Beside  the 
symmetry  plane,  also  the  symmetry  axis  shows  up  with  fast  velocity  in  the 
lower  part  of  the  figure.  The  polarization  deviations  on  the  right-hand  side 
of  figm-e  3.12b  are  given  for  0°,  10°,  ...,  50°incidence  and  may  be  compared 
with  the  data  in  figure  3.7.  The  major  features  of  the  data  set  are  explained 
by  the  model. 

The  results  for  assumed  dry  and  saturated-crack  model  are  quite  different. 
The  saturated-crack  model  gave  reasonable  fit  to  the  data  and  a  strike  angle 
of  the  high-velocity  plane  of  118°,  while  the  dry-crack  was  in  clear  conflict 
with  the  data. 

Now  we  perform  the  more  general  nonlinear  search  for  the  parameter  set 
(6,  $,  77,  r),  which  unhke  the  previous  example  does  not  assume  a  particular 
physical  fracture  model.  With  4  parameters  and  the  given  boimds  on  77  and  r 
this  search  is  not  particularly  compute-intensive,  but  guarantees  the  globally 
optimal  solution,  given  that  the  botmds  are  appropriate  and  the  sampling 
is  sufficiently  fine.  The  misfit  pattern  in  figure  3.13a  is  not  unhke  the  one 
for  the  sat-urated  crack  model.  However,  the  optimum  in  the  Southwest  is 
more  pronoimced  now  and  the  variance  reduction  is  three  percent  larger, 
now  about  44%.  The  resulting  parameters  are  0  =  203.3°,  $  =  41.05°, 
77  =  0.887  and  r  =  1.09.  The  maximum  P-velocity  difference  of  the  model 
is  about  4.7%.  The  change  in  orientation  from  the  fit  of  model  HCSl  by  a 
few  degrees  is  insignificant  considering  the  large  confidence  region  in  figure 
3.12a.  In  the  general  model  the  confidence  region  is  substantieilly  smaller. 
For  the  polarizations,  the  main  difference  is  the  behaviom  in  the  South- 
Western  part,  where  the  general  model  predicts  a  somewhat  smaller  effect. 
This  corresponds  to  the  smaller  velocity  contrast  in  that  region.  The  high- 
velocity  symetry  plane  strikes  with  about  113°  and  dips  with  about  49°  to 
the  Northeast. 

In  the  search  we  varied  77  and  r  independently.  But  can  they  be  resolved 
independently  by  the  data  constraints  or  are  they  subject  to  a  major  tradeoff? 
Figure  3.14  shows  the  misfit  around  the  optimum  model  keeping  the  symetry 
axis  orientation  fixed  and  varying  77  and  r.  The  lowest  contour  again  shows 
the  confidence  region  of  admissible  77  and  r,  which  extends  for  77  from  0.76 
to  1.07  and  for  r  from  1.  to  1.2  .  The  isotropic  case  at  77  =  1  and  r  =  1  is 
well  outside  the  highest  contour  for  four  times  the  critical  misfit  separation. 

The  strong  ascent  for  A'  near  0  strongly  discriminates  against  the  dry 
crack  model. 
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3.6  WHAT  CAUSES  THE  EFFECTIVE  ANISOTROPY? 


In  principle,  intrinsic  anisotropy  as  well  as  macroscopic  effects  from  ori¬ 
ented  cracks,  fractures  or  layering  can  cause  anisotropy  of  this  size.  To 
rmderstand  the  natme  of  this  anisotropy,  it  is  of  large  value  to  consider  ge¬ 
ological  evidence  for  the  region.  The  array  is  located  on  outcrop  of  granite 
and  gneiss.  While  the  granite  apparently  has  no  preferred  alignment,  the 
gneiss  shows  strong  foliation  with  consistent  orientations.  The  strike  is  ap¬ 
proximately  120° and  the  dip  50-60°(0tt,  pers.comm.).  This  closely  coincides 
with  the  orientation  of  the  high-velocity  planes  for  the  best-fitting  model  from 
the  imconstrained  inversion  (figure  3.13)  suggesting  that  rock  foliation  is  the 
major  cause  of  anisotropy  in  this  case.  The  size  of  the  anisotropy  suggests 
that  the  foliation  direction  is  spatially  consistent  over  some  distance,  perhaps 
several  kilometers.  In  fact,  surface  geology  does  suggest  that  this  is  indeed 
the  case  (Ott,  pers.comm.).  Rock  foliation  as  a  cause  of  effective  anisotropy 
has  been  discussed  before  (f.e.  Liischen  et  al.,  1991). 

But  is  rock  foliation  the  only  factor  involved  and  how  does  the  overall 
anisotropy  relate  to  intrinsic  properties  of  minerals?  Gneisses  typically  con¬ 
tain  large  amounts  of  mica,  a  mineral  giving  rise  to  very  large  anisotropy. 
Since  these  mineral  are  usually  quite  well-oriented  (Wang  et  al.,  1975),  we 
may  expect  major  effect  on  the  observed  anisotropy.  In  fact,  mica  orienta¬ 
tions  are  typically  dominant  in  forming  macroscopic  gneiss  foliation  planes 
(Suppe,  1985).  Pure  mica,  however,  does  not  explain  the  observed  effect. 
Aleksandrov  and  Ryzhova  (1961)  give  elastic  constants  for  muscovite,  for 
which  we  compute  77  =  0.625  and  r  =  3.88  .  The  latter  value  is  clearly 
incompatible  with  our  results.  More  realistic  models  would  include  struc¬ 
tural  complexity,  the  occurence  of  several  types  of  intrinsic  anisotropy,  and 
fractures  or  cracks  in  the  subsurface.  In  fact,  we  have  physical  models  at 
our  disposition  predicting  some  of  these  effects.  For  a  medium  consisting 
of  anisotropic  layers,  Schoenberg  and  Doiima  (1988)  give  expressions  for  av¬ 
eraging  elastic  constants  to  obtain  effective  elastic  constants  in  the  long- 
wavelength  assumption.  For  such  a  mediiun,  also  effective  parameters  fj  and 
f  are  composed  of  the  averaged  elastic  constants  as 


V 


Cl3  +  C44 
C33  —  C44 


(26) 


f  =  ^ ^  (27) 

C33  —  C44 


using 


Cii  =<  Cii  >  —  <  Ciz^ /C33  >  +  <  C13/C33  >^  /  <  I/C33  >  (28) 

ci3  =<  C13/C33  >  /  <  I/C33  >  (29) 

C33  =  1/  <  I/C33  >  (30) 
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C44  =  1/  <  1/C44  > 


(31) 


where  <>  denotes  the  thickness- weighted  average  of  the  respective  quantity. 
If  we  take  the  simple  approach  of  assuming  that  due  to  their  large  effect  only 
mica  contribute  to  the  bulk  anisotropy,  we  may  simply  compute  r}  and  r  for 
different  percentages  of  mica  content.  For  a  mica  fraction  of  30%  mica  after 
Wang  et  al.  (1975)  and  70%  other  material  assumed  to  be  isotropic  we  obtain 
77  =  0.816  and  r  =  1.74  .  While  this  is  far  outside  the  confidence  region  of 
the  observation,  values  for  3%  to  8%  mica  content  fall  within  the  confidence 
region  (figure  3.14).  These  are  low  values  compared  with  mica  contents  in 
typical  gneisses  (Wang  et  al.,  1975).  We  have  to  keep  in  mind  though  that 
certainly  not  all  of  the  mica  minerals  are  fully  aligned.  Other  constituents 
may  also  play  a  minor  role  in  defining  the  bulk  anisotropy. 

Satisfying  the  polarization  data,  this  gneiss  model  is  also  a  plausible  ex¬ 
planation  in  conjimction  with  the  geological  data.  Another  possible  model, 
however,  is  that  of  cracks  or  fractures  aligned  in  the  direction  of  planes  of 
weakness  which  we  coincide  with  the  gneiss  foHation.  Fitting  a  model  of 
aligned  cracks  after  Hudson  (1981)  had  given  essentially  the  same  orienta¬ 
tion.  Resulting  parameters  were  strike  114°  and  dip  50°  for  the  symetry 
planes.  We  also  obtained  a  crack  density  of  e  =  0.09  and  two  Lame  parame¬ 
ters  of  the  crack  material  A'  =  0.049  GPa  and  /x'  =  0.0  GPa.  The  associated 
error  bars  discriminate  clearly  against  the  dry-crack  case,  where  =  O.GPa 
.  Such  dry  cracks  give  rise  to  large  values  of  r,  which  we  have  previously 
argued,  are  in  conflict  with  the  data.  Based  on  the  achieved  misfit,  the  model 
of  cracks  opening  along  the  fofiation  planes  can  not  be  ruled  out. 

It  appears  difficult  to  distinguish  between  models  of  intrinsic  anisotropy 
and  cracks  or  fractures.  In  this  case,  one  should  consider  all  possible  models. 
In  fact,  the  occurence  of  cracks  along  planes  of  weakness  given  by  foliation 
planes  is  not  unlikely.  For  small  pressures  in  the  laboratory,  it  is  very  com¬ 
mon  to  observe  opening  of  cracks  along  the  foliation  direction.  This  can  be 
demonstrated  f.e.  with  the  data  of  Wang  et  al.  (1975).  From  their  data  we 
obtain  values  of  r  w  1.25  at  Ikbar  pressure  and  much  larger  values  up  to 
10  at  lower  pressure.  Opening  of  cracks  can  be  caused  by  the  local  stress 
field.  Values  for  the  strike  of  the  maximiun  horizontal  stress  direction  an  in 
the  area  are  around  149°  in  the  KTB  borehole  (Brudy,  Fuchs  and  Zoback, 
1993),  but  there  are  also  values  in  the  range  of  115°  (f.e.  Falkenberg  granite, 
Rummel  et  al.,  1983).  These  values  are  not  far  from  the  foliation  strike. 
Hence  crack  opening  along  the  foliation  direction  is  not  unlikely,  particularly 
if  the  fiuid  pressure  were  high.  However,  we  note  that  anisotropy  in  our  data 
set  appears  to  be  due  to  foliation  rather  than  vertical  cracks  perpendicular 
to  the  smallest  horizontal  stress  direction.  On  the  other  hand,  the  regional 
stress  field  may  have  an  effect  due  to  opening  of  cracks  preferentially  parallel 
to  the  planes  of  weakness,  which  are  given  by  the  gneiss  foliation  planes. 

Perhaps,  more  insight  can  be  gained  if  also  shear  waves  are  used.  Figure 
3.15  gives  the  prediction  of  the  shear  wave  velocities  for  the  optimum  model. 
However,  unambiguous  determination  of  upper  crustal  anisotropy  requires  in- 
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cidence  angles  within  the  shear  wave  window  and  sources  close  enough  to  the 
receivers  to  not  be  sensitive  to  more  distant  anisotropy.  These  conditions  are 
difficult  to  meet  for  our  region.  In  contrast,  for  the  study  of  local  anisotropy 
P-wave  polarization  data  may  be  used  from  local,  regional  and  teleseismic 
events,  if  the  3-component  array  meets  the  measurement  requirements. 
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TABLE  3.1:  Events  Used  in  this  Study 

(Location  and  time  from  PDE  Monthly  Bulletin,  if  not  stated  otherwise) 


Year,  Day,  Hour,  Min,  Sec. 

Latitude  Longitude  Depth 

Phase 

1990  123  01  03  38 

43.315 

19.890 

5. 

Pn 

1990  240  20  21  22 

36-267 

27.218 

39.9 

Pn 

1990  246  10  48  33 

45-915 

15.873 

20. 

Pn 

1990  253  12  17  36 

70-770 

-13.803 

10. 

P 

1990  258  23  07  43 

64-655 

-17.617 

10. 

P 

1990  285  17  30  00 

37-250 

-116.490 

0. 

P 

1990  297  15  04  14 

73-360 

54.670 

0. 

P 

1990  318  18  11  58 

-22-258 

-138.805 

0. 

PKP 

1990  325  16  59  58 

-21.9 

-138.980 

0. 

PKP 

1990  331  04  37  58 

43.853 

16.633 

24. 

Pn 

1990  331  04  51  36 

43.895 

16.641 

10. 

Pn 

1990  347  00  24  26 

37.3 

15.438 

11.1 

Pn 

1990  348  03  21  27 

39.347 

15.355 

276. 

Pn 

1990  350  15  45  41 

41.361 

43.715 

33. 

P 

1990  355  13  12  53 

-18.891 

-177.971 

457. 

PKP 

1991  038  07  12  48 

47.590 

15.503 

10. 

Pn 

1991  042  15  43  44 

44.871 

6.704 

14. 

Pn 

1991  059  15  29  42 

51.690 

16.360 

0. 

Pn 

1991  074  03  24  09 

34.343 

26.389 

7.1 

Pn 

1991  078  02  51  26 

39.260 

20.429 

10. 

Pn 

1991  094  19  00  00 

37.296 

-116.313 

0. 

P,PcP 

1991  100  01  08  40 

37.359 

36.221 

10. 

P 

1991  111  08  51  56 

37.972 

19.952 

39.8 

Pn 

1991  116  22  35  11 

43.446 

16.244 

12. 

Pn 

1991  119  09  12  48 

42.453 

43.673 

17.2 

P 

1991  119  18  30  42 

42.503 

43.899 

14.3 

P 

1991  120  03  40  36 

51.690 

16.206 

10. 

Pn 

1991  122  10  15  17 

47.929 

16.209 

10. 

Pn 

1991  123  23  41  02 

42.647 

43.263 

10.8 

P 

1991  136  02  06  17 

52.309 

7.649 

52.4 

Pn 

1991  138  17  14  59 

-21.832 

-139.014 

0. 

PKP 

1991  143  19  42  56 

51.416 

15.850 

10. 

Pn 

1991  149  18  59  58 

-22.256 

-138.794 

0. 

PKP 

1991  149  20  24  40 

45.016 

8.213 

10. 

Pn 

1991  154  10  22  40 

40.048 

42.859 

27.8 

P 

1991  165  17  59  58 

-21.944 

-138.988 

0. 

PKP 

1991  166  00  59  20 

42.461 

44.009 

9.4 

P 

1991  185  06  26  32 

42.387 

44.116 

20.4 

P 

1991  191  23  57  20 

51.424 

16.217 

0. 

Pn 

1991  196  18  09  58 

-21.877 

-138.963 

0. 

PKP 

1991  200  01  19  52 

45.344 

21.123 

10. 

Pn 
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1991  200  01  27  32 

45.312 

21.053 

10. 

Pn 

1991  222  05  23  48 

51.428 

16.242 

0. 

Pn 

1991  252  18  36  52 

51.414 

16.220 

0. 

Pn 

1991  257  19  00  00 

37.226 

-116.428 

0. 

P,PcP 

1991  279  01  46  48 

41.096 

43.409 

18.3 

P 

1991  281  03  31  16 

45.587 

149.049 

146. 

P 

1991  291  19  12  00 

37.063 

-116.045 

0. 

P,PcP 

1991  301  00  21  32 

44.265 

21.456 

67. 

Pn 

1991  303  16  24  21 

45.029 

9.937 

10.7 

Pn 

1991  304  09  31  17 

45.011 

10.059 

10. 

Pn 

1991  324  01  54  17 

46.778 

9.519 

10. 

Pn 

1991  325  02  16  32 

45.491 

21.176 

27. 

Pn 

1991  326  07  40  01 

45.019 

9.985 

11. 

Pn 

1991  327  01  06  31 

51.428 

16.243 

0. 

Pn 

1991  332  17  19  56 

36.924 

49.603 

15.7 

P 

1991  336  08  49  40 

45.498 

21.115 

9. 

Pn 

1991  336  09  04  43 

45.437 

21.247 

10. 

Pn 

1991  337  16  58  23 

44.994 

9.965 

14. 

Pn 

1991  338  13  31  31 

44.978 

9.932 

19.6 

Pn 

1991  351  06  38  17 

47.393 

151.499 

157. 

P 

1991  353  03  12  22 

45.907 

21.569 

10. 

Pn 

1991  353  18  55  17 

28.102 

57.304 

27. 

P 

1991  354  02  06  05 

24.720 

93.103 

41. 

P 

1991  356  08  43  13 

45.533 

151.021 

25. 

P 

1992  001  10  12  20 

44.984 

9.962 

23.1 

Pn 

1992  002  19  41  45 

5.638 

-73.832 

134. 

P 

1992  014  12  22  22 

51.534 

16.181 

10. 

Pn 

1992  023  06  27  39 

38.356 

20.307 

10. 

Pn 

1992  033  00  31  30 

-51.547 

139.704 

10. 

PKP 

1992  048  00  01  57 

79.191 

124.482 

10. 

P 

1992  048  19  23  14 

46.761 

9.505 

10. 

Pn 

1992  052  20  50  32 

45.463 

14.328 

11. 

Pn 

1992  062  12  29  40 

52.915 

159.886 

39. 

P 

1992  065  14  39  10 

52.9 

159.619 

45. 

P 

1992  067  01  53  38 

10.210 

-84.323 

79. 

P 

1992  073  17  18  43 

40.2 

39.500 

27.2 

P 

1992  103  08  19  30 

43.964 

13.022 

10. 

Pn 

1992  104  01  20  01 

51.157 

5.815 

15. 

Pn 

1992  112  22  28  06 

43.264 

17.961 

33. 

Pn 

1992  126  16  48  03 

43.286 

17.471 

10. 

Pn 

1992  129  06  44  39 

47.248 

9.584 

10. 

Pn 

1992  129  07  51  25 

47.233 

9.564 

10. 

Pn 

1992  130  05  37  58 

49.953 

7.409 

10. 

Pn 

1992  136  00  43  42 

47.239 

9.543 

5. 

Pn 

1992  140  14  25  04 

49.105 

6.623 

18. 

Pn 
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1992  142  04  59  57 

41.591 

88.805 

0. 

P 

1992  148  18  14  00 

44.406 

11.737 

10. 

Pn 

1992  149  09  52  42* 

44.150 

11.466 

5. 

Pn 

1992  149  12  26  02* 

43.966 

11.977 

13.6 

Pn 

1992  286  13  09  56 

29.888 

31.223 

24.7  P 

1992  297  23  19  47 

42.503 

45.073 

33. 

P 

1992  323  21  10  41 

38.300 

22.428 

10. 

Pn 

1992  326  05  07  23 

35.931 

22.443 

70.2  P 

1993  027  23  41  00 

35.963 

22.511 

62.8  Pn 

*  Location  and  time  information  from  ROM 
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b)  GERESS  Array  Configuration 

622000  623000  624000  625000  626000  627000 


Figure  3.1.  a)  Geographic  location  of  the  GERESS  array  in  Southeastern 
Germany,  b)  Array  locations  in  Gauss-Kruger  coordinates.  Data  from 
25  vertical-component  shortperiod  instruments  and  three-component 
instruments  (triangles  and  rectangles)  are  used  in  this  paper.  The 
topographic  variation  across  the  array  is  about  200  meters. 
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POLARIZATION  ATTRIBUTES  FOR  A  REGIONAL  EVENT 


t+01  X  t+Ol  X  t+OF  X  8+01  X 
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igure  d,2:  rolanzation  attributes  oi  three-component  data  indicative  tor 
body  wave  phases.  The  three-component  amplitude,  the  root-mean- 
square  amplitude  of  the  three-component  seismogram  shows  the  main 
phases  clearly.  Rectilinearity  and  component  ratios  allow  determina¬ 
tion  of  the  wave  type. 


Apparent  Incidence  Angle(deg) 


Free-Surface  Effect  on  Polarization  Inc.AngIe 


Figure  3.4:  Free-surface  effect  on  polarization  incidence  angles  for  a  set  of 
Poisson  ratios. 
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array  of  3-component  stations 

V  V  V  -V  V  V  ^ - 

predicted  propagation  direction 


Figure  3.8:  Illustration  of  plane  P-wave  propagation  in  the  crust  under  the 
three-component  receiver  array.  In  an  anisotropic  medium,  the  po¬ 
larization  directions  generally  deviate  from  the  propagation  directions. 
This  deviation,  here  shown  for  the  incidence  angle  is  used  to  infer 
effective  crustal  anisotropy.  Note  that  the  propagation  direction  typi¬ 
cally  deviates  from  the  predicted  propagation  direction,  due  to  source 
mislocation  and  distant  heterogeneity.  With  observations  of  polariza¬ 
tion  and  propagation  directions  distant  effects  can  be  excluded.  Local 
disturbances  of  the  wavefront  are  downweighted  by  averaging  over  all 
three-component  stations. 
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The  Dry  Crack  Model  (Lower  Hemisphere  Projections) 
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Figure  3.9:  Lower  hemisphere  representation  of  P-wave  velocity  and  polariza¬ 
tion  anomalies  for  (tj,  r)  =  (1-43,  2.13)  corresponding  to  the  cry-crack 
model  HCDl  of  Crampin  (1984). 
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4.453 


The  Saturated  Crack  Model  (Lower  Hemisphere  Projections) 


Fimre  3  10-  Lower  hemisphere  representation  o{  P-wave  velocity  and  po- 
larisation  anomaUes  for  (,.  r)  =  (0.82  LOO)  corresponding  to 
saturated-crack  mode!  HCSl  of  Crampin  (1984). 
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^igure  3-11:  Fitting  of  tiie  dxy-crack  model.  Misnt  on  tiie  spFere  (23)  is 
shown  in  degree^,  depending  on  the  sjmxaetrj'  axis  orientation.  Smaller 
misnts  axe  found  near  for  axes  located  near  the  horizontal  plane.  Tnis 
drv-crack  model  doesn't  improve  the  misnt  and  hence  can't  explain  the 


Misfit  (lopondinj*  on  Syin 

'I'lU'lii: 'iOK.O  IMii  ^36.0  I'.la:  I 


Figure  3.12:  Fititiug  of  the  sat'orated-crack  model:  a)  misni  in  degree- 
depending  on  symetry-  axis  orientation  from  O^to  90-:ncidence.  b) 
best-fitting  model  shown  in  a  partial  lower  hemispiiere  display"  (0®to 
55° incidence).  The  symmetry  axis  has  strike  208°  and  dip  36°  to  the 
South.  The  lowest  contotir  line  giv'es  the  critical  level.  Polarization 
deviations  are  given  for  0°,  10°,—,  50°incidence. 
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Misi'i(  (l(‘|)(!!uliii}»  on  Symelry  Axis  Oriciilation 

rholii: 1'hi.:.|I.OS  K(ii::().«K7TiUi::l.00S-.:.1H.K202 
(Idwcr  lieinisphcrc  pnijoction) 


Misfit  depending  on  Eta  and  Tau 
Theta=205300  Phi=4 1.0500  et2=0.SS7000  taa=L09000  S=38.S202 


i  i 

'02.-:  172 

91.5755 
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49.4085 

38.8668 


Figure  3.14;  MisSt  around  the  optimxmi  model  (black  square)  keeping  the 
orientation  fixed.  The  lo-west  contour  line  shows  the  confidence  region 
around  the  best  fitting  model.  Numbers  0  to  13  refer  to  percentages  o 
mica  content  in  a  gneiss  model  (see  text). 
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Shear-Wave  Velocities  for  the  Optimum  Model 
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Figure  3,15:  Prediction  of  shear  wave  velocities  for  the  optirnuin  model  on 
the  full  lower  hemisphere  (O'* to  90®).  qSP  is  polarised  paralkd  to  the 
symetry  plane  and  qSR  perpendicular. 
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